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1 INTRODUCTION 



1.1 Overview 

The field of neutron physics has become an integral part of investigations into an array of important issues 
that span fields as diverse as nuclear and particle physics, fundamental symmetries, astrophysics and cos- 
mology, fundamental constants, gravitation, and the interpretation of quantum mechanics. The experiments 
employ a diversity of measurement strategies and techniques - from condensed matter and low temperature 
physics, optics, and atomic physics as well as nuclear and particle physics - and address a wide range of 
issues. Nevertheless, the field possesses a certain coherence that derives from the unique properties of the 
neutron as an electrically neutral, strongly interacting, long-lived unstable particle that can be used either 
as the probe or as an object of study. This review covers some of the important new contributions that 
neutrons have made in these diverse areas of science. By "fundamental" neutron physics we mean that class 
of experiments using slow neutrons that primarily address scientific issues associated with the Standard 
Model (SM) of the strong, weak, electromagnetic, and gravitational interactions and their connection with 
issues in astrophysics and cosmology. 

Neutrons experience all known forces in strengths that make them accessible to experimentation. It is 
an amusing fact that the magnitude of the average neutron interaction energy in matter, in approximately 
1 T magnetic fields, and over a one meter fall near the surface of the Earth is the same order of magni- 
tude (« 100 neV). This coincidence leads to interesting and occasionally bizarre strategies for experiments 
that search for gravitational effects on an elementary particle. The experiments include measurement of 
neutron-decay parameters, the use of parity violation to isolate the weak interaction between nucleons, and 
searches for a source of time reversal violation beyond the SM. These experiments provide information that 
is complementary to that available from existing accelerator-based nuclear physics facilities and high-energy 
accelerators. Neutron physics measurements also address questions in astrophysics and cosmology. The the- 
ory of Big Bang Nucleosynthesis needs the neutron lifetime and the vector and axial vector weak couplings 
as input, and neutron cross sections on nuclei are necessary for a quantitative understanding of clement 
creation in the universe. 

Free neutrons are unstable with a 15 minute lifetime but are prevented from decaying while bound in 
nuclei through the combined effects of energy conservation and Fermi statistics. They must be liberated from 
nuclei using nuclear reactions with MeV-scale energies to be used and studied. We define "slow" neutrons to 
be neutrons whose energy has been lowered well below this scale. The available dynamic range of neutron 
energies for use in laboratory research is quite remarkable, as shown in Table [T] Thermodynamic language 
is used to describe different regimes; a neutron in thermal equilibrium at 300 K has a kinetic energy of only 
25meV. Since its de Broglie wavelength (0.18 nm) is comparable to interatomic distances, this energy also 
represents the boundary below which coherent interactions of neutrons with matter become particularly 
important. The most intense sources of neutrons for experiments at thermal energies are nuclear reactors, 
while accelerators can also produce higher energy neutrons by spallation. 

Neutron decay is an important process for the investigation of the Standard Model of electroweak in- 
teractions. As the prototypical beta decay, it is sensitive to certain SM extensions in the charged-current 
electroweak sector. Neutron decay can determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\V u d\ through increasingly precise measurements of the neutron lifetime and the decay correlation coefficients. 

Searches for violations of time-reversal symmetry and/or CP symmetry address issues which lie at the 
heart of cosmology and particle physics. Among the important issues that can be addressed by neutron 
experiments is the question of what mechanisms might have led to the observed baryon asymmetry of 
the universe. Big Bang cosmology and the observed baryon asymmetry of the universe appear to require 
significantly more T-violation among quarks in the first generation than is predicted by the Standard Model. 
The next generation of neutron electric dipole moment (EDM) searches, which plan to achieve sensitivities 
of 10 _27 e • cm to 10~ 28 e • cm, is the most important of a class of experiments aiming to search for new 
physics in the T- violating sector. 

The last decade has also seen qualitative advances in both the quantitative understanding of nuclei, 
especially few body systems, and in the connection between nuclear physics and quantum chromodynamics 
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Term Energy Velocity Wavelength Temperature 

(m/s) (nm) (K) 



ultracold < 0.2 ^eV 

very cold 0.2 [icV < E < 50 /xeV 

cold 0.05 mcV < E < 25 mcV 
thermal 25 meV 

epithermal 25 meV < E < 500 keV 
fast > 500keV 



< 6 
6 < v < 100 
100 < v < 2200 
2200 

2200 < v < 1 x 10 7 
> 1 x 10 7 



> 64 
4 < A < 64 
0.18 < A < 4 
0.18 



< 0.002 
0.002 < T < 0.6 
0.6 < T < 300 
300 



Table 1: Common terminology and spectrum of neutron energies. 



(QCD). Low energy properties of nucleons and nuclei, such as weak interactions in n-A systems, low energy 
n-A scattering amplitudes, and the internal electromagnetic structure of the neutron (its electric polariz- 
ability and charge radius) are becoming calculable from the SM despite the strongly interacting nature of 
these systems. These theoretical developments are motivating renewed experimental activity to measure un- 
determined low energy properties such as the weak interaction amplitudes between nucleons and to improve 
the precision of other low energy neutron measurements. The ultimate goal is to illuminate the strongly 
interacting ground state of QCD, the most poorly-understood sector of the SM, and to clarify the connection 
between the strong interaction at the hadron level and the physics of nuclei. 

This review presents and discusses the status of the experimental efforts to confront these physics ques- 
tions using slow neutrons. The improvements in precision required to address these questions are technically 
feasible and have spurred both new experimental efforts and the development of new neutron sources. We 
also discuss some of the new proposed facilities under construction. It is not possible to cover the large 
volume of work in a review of this scope, so we refer the reader to a number of more specialized reviews 
wherever appropriate. Instead, we emphasize recent experiments and those planned for the near future. 
Neutron experiments form part of a larger subclass of low energy precision measurements which test the 
SM [1-5]. There are texts that cover a broader survey of topics and provide historical context [6-9]. Tests 
of quantum mechanics using neutron intcrfcrometry are not discussed here but are covered in detail in a 
recent comprehensive text [10]. Proceedings of recent international conferences also describe many planned 
projects or measurements in progress [11,12]. 

1.2 Neutron Optics 

A basic introduction to neutron optics is very useful for understanding how experiments using cold and 
ultracold neutrons are conducted. We will use it as a unifying concept in this review, since a large number of 
measurements possess a transparent neutron optical interpretation. For a rigorous development of neutron 
optics theory, see Sears [13]. 

An optical description for particle motion in a medium is possible for that subset of amplitudes in which 
the state of the medium is unchanged. For these events the wave function of the particle is determined by a 
one-body Schrodingcr equation 



where ip(r) is the coherent wave and v(r) is the optical potential of the medium. The coherent wave and 
the optical potential satisfy the usual coupled system of equations (Lippmann-Schwingcr) of nonrelativistic 
scattering theory. To an excellent approximation the optical potential for slow neutrons is local and energy- 
independent and given by 



2m 



A + v(r) tp(r) = Eip(r), 



(1) 




(2) 
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where Ni is the number density of scatterers and bi is the coherent scattering length for element I. All 
neutron-matter interactions contribute to the scattering length: the dominant contributions come from 
the neutron-nucleus strong interaction and the interaction of the neutron magnetic moment with magnetic 
fields. Both interactions are spin-dependent and therefore the optical potential possesses two channels with 
J = J± 1/2. The coherent scattering length b is the sum of the scattering lengths in both scattering channels 
weighted by the number of spin states in each channel. From a quantum mechanical point of view, this is the 
total amplitude for a neutron to scatter without a change in the internal state of the target. The effect of 
the optical potential for a nonabsorbing uniform medium with positive coherent scattering length is to slow 
down the neutrons as they encounter the potential step due to the matter, thereby decreasing the neutron 
wave vector, K, within the medium. A neutron index of refraction can be defined by this relative change in 
the magnitude of the wavevector n = K/k. Conservation of energy at the boundary determines the relation 
to the optical potential 

n 2 = 1 - v /E. (3) 

Incoherent effects due to neutron absorption and incoherent scattering with the medium remove probability 
density from the coherent wave and can be modeled by adding an imaginary part to the optical potential 
whose magnitude is determined by the optical theorem of nonrelativistic scattering theory. The scattering 
lengths therefore become complex. 

It may seem counterintuitive that the attractive neutron-nucleus interaction can give rise to a repulsive 
neutron optical potential. For a weak potential there is indeed a one-to-one relation between the sign of the 
zero-energy scattering phase shift and the sign of the potential, but the neutron-nucleus interaction is strong 
enough to form bound states. If one recalls Levinson's theorem from nonrelativistic scattering theory, which 
says that the zero energy phase shift increases by tt for each bound state, one can see that it is possible for an 
attractive potential to produce a negative phase shift. Simple models [14] and more detailed treatments [15] 
of the neutron-nucleus interaction show that for most nuclei the scattering lengths are indeed positive. A 
summary of all neutron scattering lengths up to 1991 [16] and a set of recommended values [17] exist in the 
literature. 

The fact that the neutron optical potential is repulsive for most materials has a number of far-reaching 
consequences for slow neutron technology that are briefly mentioned here. Since the neutron index of 
refraction is less than unity, the well-known light optical phenomenon of total internal reflection for rays 
at an interface becomes total external reflection for neutrons. Neutron guides [18], the neutron analog of 
optical fibers, are used to conduct neutron beams far from the neutron source by total external reflection from 
mirrors with negligible loss in intensity within their phase-space acceptance. The phase-space acceptance 
can be increased over that of a uniform medium by multilayer coatings called supermirrors, which expand 
the effective critical angle for total external reflection through constructive interference of scattering from 
different layers. 

If the mirror is made of magnetic material with a magnetic scattering amplitude of the same size as the 
nuclear scattering amplitude, the optical potential can be large for one neutron spin state and close to zero 
for the opposite spin state. In this case the reflected neutron beam is highly polarized. Supermirror neutron 
polarizers which polarize essentially all of the phase space of a cold neutron beam transported by a neutron 
guide to greater than 98 % are routinely available. A polarized neutron beam can also be produced through 
transmission of an unpolarized neutron beam through a polarized medium; both polarized hydrogen (using 
spin-dependent scattering) and polarized 3 He (using spin-dependent absorption) have been used. 

Neutrons whose kinetic energy is lower than the optical potential of the medium that are incident on 
a surface will be reflected at all angles of incidence. Such neutrons can therefore be trapped in a material 
bottle through total external reflection and are called ultracold neutrons (UCN). The possibility of performing 
measurements on an ensemble of free neutrons almost at rest in the lab is a decisive advantage in certain 
experiments. Neutrons in this energy range can also be confined by magnetic field gradients that act on the 
neutron magnetic moment. In addition, such neutrons can only rise a meter or so in the gravitational field 
of the Earth. In both of these cases the "optical potential" is simply the normal potential energy of the 
neutron in the presence of the relevant external field. 

The expression for the low energy neutron scattering length of an atom away from nuclear resonances for 
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unpolarized atoms and neutrons and nonmagnetic materials is 

b = b nuc + Z(b ne + b s )[l - f(q)] + b pol , (4) 

where b nuc is the scattering amplitude due to the neutron-nucleus strong force, b ne is the neutron-electron 
scattering amplitude due to the internal charge distribution of the neutron, b s is the Mott-Schwinger scatter- 
ing due to the interaction of the magnetic moment of the neutron with the v x E magnetic field seen in the 
neutron rest frame from electric fields, b po i is the scattering amplitude due to the electric polarizability of the 
neutron in the intense electric field of the nucleus, and f(q) is the charge form factor (the Fourier transform 
of the electric charge distribution of the atom) . The electromagnetic contribution to the scattering lengths 
from b ne and b s are exactly zero for forward scattering due to the neutrality of the atoms, which forces the 
charge form factor f(q) — > 1 as q — > [13]. The weak interaction also makes a contribution to the scattering 
length which possess a parity-odd component proportional to s* n ■ p, where is the neutron spin and p is 
the neutron momentum. The relative sizes of these contributions to the scattering length from the strong, 
electromagnetic, and weak interactions are roughly 1 : 10~ 3 : 10~ 7 . 

1.3 Neutron Sources 

Most fundamental neutron physics experiments are conducted with slow neutrons for two main reasons. 
First, slower neutrons spend more time in an apparatus. Second, slower neutrons can be more effectively 
manipulated through coherent interactions with matter and external fields. Free neutrons are usually created 
through either fission reactions in a nuclear reactor or through spallation in high Z targets struck by GeV 
proton beams. We briefly examine these neutron sources and the process by which cold and ultracold 
neutrons (UCN) are produced starting from neutrons with energies several orders of magnitude greater. 

Neutrons are produced from fission in a research reactor at an average energy of 2 MeV. They are slowed 
to thermal energy in a moderator (such as heavy or light water, graphite, or beryllium) surrounding the 
fuel. The peak core fluence rate of research reactors is typically in the range 10 14 cm~ 2 s _1 to 10 15 cm~ 2 s _1 . 
To maximize the neutron density, it is necessary to increase the fission rate per unit volume, but the power 
density is ultimately limited by heat transfer and material properties. In the spallation process, protons 
(typically) are accelerated to energies in the GeV range and strike a high Z target, producing approximately 
20 neutrons per proton with energies in the fast and epithermal region [19]. This is an order of magnitude 
more neutrons per nuclear reaction than from fission. Existing spallation sources yield peak neutron rates 
of lO^s -1 and 10 17 s _1 . Although the time -averaged fluence from spallation neutron sources is presently 
about an order of magnitude lower than for fission reactors, there is potentially more room for technical 
improvements in the near-term future due in part to the relaxation of the constraints needed to maintain a 
nuclear chain reaction [20] and to the potential for future developments in GeV proton accelerator technology. 

The main feature that differentiates spallation sources from reactors is their convenient operation in 
a pulsed mode. At most reactors one obtains continuous beams with a thcrmalized Maxwellian energy 
spectrum. In a pulsed spallation source, neutrons arrive at the experiment while the production source is 
off, and the frequency of the pulsed source can be chosen so that slow neutron energies can be determined by 
time-of-flight methods. The lower radiation background and convenient neutron energy information can be 
advantageous for certain experiments. The frequency is typically chosen to lie in the 10 Hz to 60 Hz range, 
so that the subsequent neutron pulses from the moderator do not overlap for typical neutron spectra and 
distances to the experimental area. 

Fast neutrons reach the thermal regime most efficiently through a logarithmic energy cascade of roughly 
20 to 30 collisions with matter rich in hydrogen or deuterium. Neutrons are further cooled by a cryogenic 
neutron moderator adjacent to the reactor core or spallation target held at a temperature of w 20 K. One 
generally wants the moderator as cold as possible to increase the phase space density of the neutrons. As 
the neutron wavelengths become large compared to the atomic spacings, the total scattering cross sections 
in matter are dominated by elastic or quasielastic processes, and it becomes more difficult for the neutrons 
to thermalize. The development of new types of cold neutron moderators is therefore an important area for 
research. 
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Table 2: Some operating parameters for major cold neutron reactor-based user facilities with active funda- 
mental physics programs. Approximate fluence rates are given as neutron capture fLuence. 



Parameter 


TT T 
ILL 


TT T 
ILL 


IN lb 1 


T?T>A/T TT 




PF1 


PF2 


NG-6 


Mephisto 


Power (MW) 


58 


58 


20 


20 


Guide length (m) 


60 


74 


G8 


30 


Guide radius (m) 


4000 


4000 


CO 


460 


Guide type (m =) 


1.2 


2 


1.2 


3 


Cross section (cm 2 ) 


6 x 12 


6 x 20 


6 x 15 


5 x 11.6 


Fluence rate (xlO 9 cm _2 s _1 ) 


4 


14 


2 


7 



It is not practical to describe specific neutron facilities in any detail in this review, but we note a few where 
the bulk of research efforts have been carried out. For thirty years the most active facility for fundamental 
neutron research has been the Institut Max von Laue - Paul Langevin (ILL) in Grenoble, France [21]. Its 
58 MW reactor is the focal point of neutron beta-decay and UCN physics in the world, and its cold neutron 
beam H113 for fundamental physics, the first large ballistic supermirror neutron guide, possesses the highest 
cold neutron intensity in the world [22,23]. The new FRM-II reactor has come online in Munich, and the 
experimental program in fundamental neutron physics has commenced [24]. The most active institutions in 
the US are the National Institute of Standards and Technology (NIST) [25], which operates the NG-6 end 
position, two monochromatic neutrons beams at 0.49 nm and 0.9 nm, and a neutron interferometer facility 
on the NG-7 guide, and Los Alamos National Laboratory (LANL), which operates the FP12 beamline [26]. 
In Russia, there are significant efforts at the Petersburg Nuclear Physics Institute (PNPI), Joint Institute 
for Nuclear Research (JINR) in Dubna, and the Kurchatov Institute and in Japan at KEK. Many smaller 
sources play an essential role in the development of experimental ideas and techniques. 

In addition to the existing sources, the last decade has seen tremendous growth in the construction of 
facilities and beamlines devoted to fundamental neutron physics. Many of these new facilities are at spallation 
sources. The Paul Scherrer Institut (PSI) operates the Swiss Spallation Neutron Source SINQ [27, 28] 
as a continuous spallation source and has constructed a cold neutron beamline FUNSPIN dedicated to 
fundamental physics [29,30]. In the US, the 2MW Spallation Neutron Source (SNS) is under construction, 
and the fundamental physics beamline (FNPB) should be operational some time in 2008 [31]. The Japanese 
Spallation Neutron Source (JSNS) is in the construction phase and is also anticipated to become operational 
in 2008. A fundamental neutron physics beamline at the JSNS with a direct line-of-sight to the neutron 
moderator would be complementary to the (curved) SNS beam and would allow for experiments using pulsed 
epithermal neutrons. Tables [2] and [3] give a few of the measured (or projected) cold-neutron beam properties 
for some of the facilities with active fundamental physics programs. As a word of caution, it is not always 
straightforward to compare the cold neutron beam properties of different sources, and experimentalists must 
consider all the parameters before determining the optimal source. 

In the neutron energy spectrum from a cold moderator there is a very small fraction whose energies 
lie below the w 100 neV neutron optical potential of matter. Such neutrons are called ultracold neutrons, 
and they can be trapped by total external reflection from material media. The existence of such neutrons 
was established experimentally in the late 1960s [32,33]. The UCN facility at the ILL employs a turbine 
to mechanically convert higher energy neutrons to UCN [34]. Although the density of neutrons is bounded 
by the original phase space in the source (Liouville's theorem), this technique produces enough UCN to 
conduct a number of unique and fundamental experiments described in part below. Over the last decade 
new types of UCN converters have been developed that can increase the phase space density through the 
use of "superthermal" techniques [35] . They involve energy dissipation in the moderating medium (through 
phonon or magnon creation) and can therefore increase the phase space density of the UCN since Liouville's 
theorem docs not apply. Superfluid helium [36] and solid deuterium [37] have been used most successfully as 
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Table 3: Some operating parameters for major cold neutron spallation-source user facilities with active (and 
proposed) fundamental physics programs. Approximate time-averaged fluence rates are given as neutron 
capture fluence. Note that for some experiments the peak fluence rather than the time-averaged fluence is 
the more important parameter. 



Parameter 


SINQ 
FunSpin 


LANSCE 
FP12 


(SNS) 
(FNPB) 


(JSNS) 


Time-averaged current (mA) 


1.2 


0.1 


(1.4) 


(0.3) 


Source rep. rate (Hz) 


dc 


20 


(60) 


(25) 


Guide length (m) 


7 


8 


(15) 


(10 to 20) 


Guide radius (m) 


oo 


oo 


(117) 


(oo) 


Guide type (m =) 


3 


3 


(3.5) 


(3) 


Cross section (cm 2 ) 


4 x 15 


9.5 x 9.5 


(10 x 12) 


(10 x 10) 


Fluence rate (xlO 8 cm~ 2 s _1 ) 


8 




' (10) ' 


(5) 



superthermal UCN sources, and other possible UCN moderating media such as solid oxygen are also being 
studied [38]. The lack of neutron absorption in 4 He along with its other unique properties makes possible 
experiments in which the measurement is conducted within the moderating medium. These developments 
have led to new construction projects for UCN facilities at LANL [37], PSI [39], FRM [40], KEK [41], NC 
State [42], Mainz, and other sources. Extensive treatments of UCN physics are found in References [43,44]. 

2 NEUTRON DECAY AND STANDARD MODEL TESTS 

There are a several reviews that discuss weak interaction physics using slow neutrons in greater detail or 
provide additional information [45-47] . A recent publication emphasizes the issue of CKM unitarity [48] . A 
comprehensive review of measurements in neutron and nuclear beta decay to test the Standard Model in the 
semileptonic sector and its possible extensions along with a comparison with other probes of similar physics 
has appeared [5]. 

2.1 Theoretical Framework 

The neutron is composed of two down quarks and an up quark, and it is stable under the strong and 
electromagnetic interactions, which conserve quark flavor. The weak interaction can convert a down quark 
into an up quark through the emission of the W gauge boson. The mass difference of the neutron and 
proton is so small that the only possible decay products of the W are an electron and antineutrino with 
the release of energy distributed among all the decay products: n — > p + e~ + D e + 0.783 MeV. There are 
two other decay channels available to the neutron. There is the radiative decay mode with a photon in the 
final state and decay to a hydrogen atom and an antineutrino [49-51]. Very recently, the radiative decay 
mode has been definitively observed and the branching ratio measured at the 10% level [52]. The decay 
into a hydrogen atom has not yet been seen although a search for n — > H + v e is under investigation [53] . 
Beta-decay experiments can test the assumptions of the SM by performing precision measurements on the 
decay product energies, momenta, and their correlations with the neutron spin. 

Free neutron decay in the Standard Model is described in lowest order by a mixed vector/axial- vector 
current characterized by two coupling strengths, gy and gA, the vector and axial- vector coupling coefficients. 
Since the momentum transfers involved in neutron beta decay are small compared to the W and Z masses, 
one can write an effective Lagrangian that describes neutron decay in the Standard Model as a four-fermion 
interaction 

C int = ?^(V ll -\A ll )(v»-a»), (5) 



6 



where V M = VV7^n 



V'eT^Vv, A^ = ippj^^ipn, and a M = V'eT^TsVv are the hadronic and leptonic 



vector and axial vector currents constructed from the neutron, proton, electron, and neutrino fermion fields, 
Gf is the Fermi decay constant, V u( j is a CKM matrix element, and A is the ratio of the axial vector and 
vector couplings. 

The V-A structure for the weak currents is incorporated directly into the standard electroweak theory 
by restricting the weak interaction to operate only on the left-handed components of the quark and lepton 
fields. A more fundamental understanding of the reason for this parity-odd structure of the weak interaction 
is still lacking. There is also no physical understanding for the measured values of the CKM mixing matrix 
elements between the quark mass eigenstates and their weak interaction eigenstates. The fact that the matrix 
is unitary is ultimately a consequence of the universality of the weak interaction gauge theory. Extensions 
to the SM which either introduce non V-A weak currents or generate violations of universality can therefore 
be tested through precision measurements in beta decay. A recent reanalysis of the constraints on non 
V-A charged currents showed that improved neutron decay measurements have set new direct limits on 
such couplings, which are typically constrained at the 5% level [5]. Complementary constraints on non 
V-A charged currents in neutron beta decay from a combination of neutrino mass limits and cosmological 
arguments have recently appeared [54] 

The probability distribution for beta decay in terms of the neutron spin and the energies and momenta 
of the decay products [55] can be written 



dW cx (g 2 v + ig 2 A )F{E e ) 



E e E y E e \ E e Ed E e E v 



(6) 



where one defines 
2ir 3 h 7 



m\<* f{l + 5 B ){gv 2 + i9A 2 ) 



9A 



J<t> 



gv 

1-|A| 2 
1 + 3|A| 2 

6 = 

A _ _ JA| 2 + |A|cos< 



B = 2 
D = 2 



1 + 3|A| 2 
|A| 2 - |A| cos0 

1 + 3|A| 2 
|A| sin 6 



1 + 3 



= (885.7 ± 0.8) s neutron lifetime 

= —1.2695 ± 0.0029 coupling constant ratio 

= —0.103 ±0.004 electron-antineutrino asymmetry 

= Fierz interference 

= — 0.1173 ± 0.0013 spin-electron asymmetry 

= 0.981 ± 0.004 spin-antineutrino asymmetry 

= (-4 ± 6) x 10~ 4 T-odd triple-product. 



In these equations, which neglect small contributions such as weak magnetism and radiative corrections, 
F(E e ) is the electron energy spectrum, p e , p„, E e , and E v are the momenta and kinetic energies of the decay 
electron and antineutrino, a n is the initial spin of the decaying neutron, and cf> is the phase angle between the 
weak coupling constants qa and gv- / is the statistical rate function and 8r is a radiative correction evaluated 
for the neutron; the product /(l + Sr) = 1.71489 ± 0.00002 has been calculated with high precision [56]. 
The spin-proton asymmetry correlation coefficient C is proportional to the quantity A + B. The values 
represent the world averages as compiled by the Particle Data Group (PDG) [57]. The parameter A can 
be extracted from measurement of either a, A, or B. If the neutron lifetime r„ is also measured, gv and 
<7j4 can be determined uniquely under the assumption that D = 0. Radiative corrections to the correlation 
coefficients exist and have been evaluated recently using effective field theory (EFT) based methods [58-63] . 
Figure [1] shows the recent history of measured values of the lifetime and correlation coefficients as used by 
the PDG [57]. 

Accurate determination of the parameters that describe neutron decay can provide important information 
regarding the completeness of the three-family picture of the SM through a test of the unitarity of the CKM 
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Figure 1: A summary plot of the measurements of the neutron lifetime and correlation coefficients that are 
used in the 2006 compilation of the PDG. The shaded bands represent the ±1 a region. 



matrix, or |Vud| 2 + |Vu S | 2 + |Kib| 2 = 1< One strong motivation for more accurate measurements of neutron 
decay parameters is to measure |Vud|, the largest term in the unitarity sum. The matrix elements |V^ S | and 
|tub| are obtained from high energy accelerator experiments. 

The conserved- vector-current hypothesis implies that gv = Gp|Kid|(l + &r), where Aj( is the nucleus- 
independent radiative correction. One can determine |V^d| by measuring the Tt values of superallowed 
+ — > + (3 transitions between isobaric analog states [64] . Recently, a new method for calculating hadronic 
corrections to electroweak radiative corrections that combines high-order perturbative QCD calculations with 
a large N-based extrapolation has been applied to superallowed nuclear beta decay [65]. This calculation 
reduced the theoretical loop uncertainty by approximately a factor of two and now yields a value of |Vud| = 
0.97377 ± 0.00027 with the uncertainty dominated by theoretical nucleus-dependent radiative corrections. 
Pion beta decay is theoretically the cleanest system in which to measure [V^dl; but the small branching ratio 
has so far precluded a measurement with enough sensitivity to compete with superallowed beta decay and 
neutron decay. The latest measurement from pion beta decay gives \V u a\ = 0.9728 ± 0.0030 [66]. 

Using values of \V ns \ and |Kib| taken from the current recommendations of the PDG, one obtains |Kid| 2 + 
\V ns \ 2 + |Vub| 2 = 0.99919 ± 0.00108, which removes the more than two standard deviation discrepancy with 
unitarity [67]. The primary reason for the significant shift is a new value of |Vus| that is about 2.5% larger 
than the previous PDG evaluation. The result is based on recent experiments in semileptonic kaon decays, 
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including measurements of the branching ratios, form factors, and lifetime [68]. (There are also renewed 
theoretical investigations to extract |U US | from hyperon decay [69] and precision experiments in hadronic tau 
decays [70]. ) It still uses the value of /+(0) calculated by Leutwyler and Roos [71], but one should note 
that other calculations exist that produce numbers that differ by as much as 2%. The value of \V u b\ has also 
changed, but it has a negligible contribution to the unitarity sum. 

Neutron beta decay offers a somewhat cleaner theoretical environment for extracting gy than the super- 
allowed + — > + (3 transitions due to the absence of other nucleons (although some radiative corrections are 
common to both systems). In neutron decay the value of |Kid| is obtained from the expression 



, 1 K/\n2 



|Kd| /(1 + <5 r )t„G f (1 + AX)(1 + 3A 2 )' (?) 

where K is a constant. Marciano and Sirlin recently refined calculations of the radiative corrections [65] to 
yield 

a (4908.7 ± 1.9) a 
|Kd| = t„(1 + 3A2) • ( g ) 

Thus, with the lifetime and the parameter A, most precisely obtained from measurements of the neutron 
spin-beta asymmetry angular correlation, |Kid| can be determined uniquely. Using the PDG values of r n 
and A, the same unitarity test gives Y^i l^ui| 2 = 0.9971 ± 0.0039, consistent with unity but less precise. 
The present situation regarding unitarity is summarized in Figure [2j Although the values from the PDG 
evaluation imply good agreement among the neutron measurements, the nuclear systems, and unitarity, the 
plot also shows unresolved disagreements in the neutron lifetime value and the spin-electron asymmetry. 
Until those issues are resolved, the results from the neutron system remain in question. 

A precision determination of \V u d\ should be seen in the context of the overall effort to determine with 
high precision all the parameters of the CKM matrix. At CLEO-c it should be possible to measure the CKM 
matrix element \V c d\ to 1 % accuracy if lattice gauge theory calculations of the required form factors can 
match the expected precision of the data [74,75]. Recent developments in lattice gauge theory have produced 
a qualitative improvement in the accuracy with which certain observables can be calculated; for example, the 
latest result for V us = 0.2219 ± 0.0026 from the MILC collaboration rivals the experimental precision [76]. 
Successful lattice gauge theory calculations of form factors for charmed quark decays would make possible 
another independent check of CKM unitarity using the first column, |Vud| 2 + \ V c a\ 2 + |Vtd| 2 = !• 



2.2 Neutron Lifetime 

2.2.1 Recent Lifetime Experiments 

Two distinct strategies have been employed to measure the neutron lifetime. The first measurements were 
performed using thermal or cold beams of neutrons and measured simultaneously both the average number 
of neutrons N and the rate of neutron decays dN/dt from a well-defined fiducial volume in the beam. 
All modern experiments in neutron decay use low energy neutrons because more particles decay within a 
given volume, thus increasing the statistical power. The beam technique requires absolute knowledge of 
the efficiencies to extract the lifetime from dN/dt = — N/r n . A significant improvement in the precision 
was obtained by implementing a segmented proton trap to better define the fiducial volume [77]. The two 
more precise experiments used such a proton trap as well as a neutron detector with an efficiency that was 
inversely proportional to the neutron velocity [78,79]. These experiments are currently limited by systematic 
effects related to neutron counting, but reductions in the uncertainties are quite feasible, and an effort is 
underway to improve the precision of one of the beam experiments. 

With the advent of UCN production techniques, a second approach to lifetime measurements was devel- 
oped that confines neutrons in material "bottles" or magnetic fields and measures the number of neutrons 
remaining as a function of time. At a time t the number of neutrons is given by N(t) = 7V(0)e -t / rn , so 
by measuring N(t) at different times, one can determine the neutron lifetime. With this technique one 
avoids the necessity of knowing the absolute neutron density and detector efficiencies but other systematic 
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Figure 2: The weak coupling constants gv and qa determined from neutron decay parameters, superallowed 
+ — * + nuclear decay transitions, and CKM unitarity. The bands represent ±1 standard deviation uncer- 
tainties. For the neutron decay parameters, the uncertainties are dominated by experimental uncertainties. 
For the superallowed + — ► + nuclear decays, the uncertainty is dominated by radiative correction cal- 
culations. In addition to the recommended values of the PDG 2006 evaluation shown for neutron decay, 
the PERKEO result [72] and the lifetime measurement of Serebrov et al. [73] are shown independently for 
comparison. 



effects arise. The measured lifetime is given by t^ \ — t^ 1 + jC and includes losses from the confinement 
cavity as well as neutron decay. To isolate r/ oss , which is typically dominated by nonspecular processes in 
the neutron interaction with the walls of the confinement vessel, one measures the lifetime in bottles with 
different surface-to-volume ratios and performs an extrapolation to an infinite volume (or equivalently zero 
collision rate). These losses depend on the UCN energy spectrum, which can change during the storage 
interval, so much work has been done to understand the spectrum evolution and loss mechanisms and to find 
surface materials with lower loss probabilities. To address losses experimentally, one group simultaneously 
measured the UCN storage time and the inelastically scattered neutrons [80], thus monitoring the primary 
loss process. Because neutrons can also occupy unbound metastable orbits in the trap, it is necessary to 
"clean" the phase space of the UCN to remove them so that the remaining losses are dominated by wall 
interactions. 

Accurate measurements using each of these independent methods are important for establishing the 
reliability of the result for r n , particularly given a recent measurement that is significantly discrepant with 
the world average from the PDG 2006 evaluation. Figure [3] shows a summary of measurements over the 
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last 20 years with competitive uncertainties. (A more historical and detailed review of neutron lifetime 
experiments is found in Ref. [81].) Seven of the experiments [78,80,82-86] contribute to a current neutron 
lifetime world average of r„ = (885.7 ± 0.8) s [57]. The agreement among the measurements is very good. 
The experiments using beams do not have as much statistical influence as those using UCN, but they are 
also in good agreement. 

Also shown on Figure [3] is the most recent experiment of Serebrov et al. that produced a value of 
r n = (878.5 ± 0.8) s [73]. The result is more than six standard deviations away from the PDG average. The 
experiment used ultracold neutrons confined both gravitationally and by a material bottle and measured 
the total storage lifetime. The loss rate on the walls is given by rZ ss = i]^, where 77 is a wall-loss coefficient 
related to the probably of losing a neutron per interaction with the wall and 7 is the loss-weighted collision 
frequency. Other loss mechanisms were considered negligible. The neutron lifetime was determined by 
extrapolating to zero collision frequency. The experiment used low temperature Fomblin oil to significantly 
reduce the value of // in comparison with previous experiments. The improvement reduced the size of the 
extrapolation necessary to extract r„, thus making the experiment less sensitive to some systematic effects. 
They also used two storage vessels with different surface-to- volume ratios and were able to change the energy 
distribution of the confined UCN, both of which change the collision frequency. The group intends to make 
additional measurements with a variable-volume trap to change the collision frequency while maintaining 
the same trap surface and vacuum conditions. 

The substantial difference between the neutron lifetime of PDG average and that of Serebrov et al. is 
not understood. It is essential to resolve the disagreement, and it can only be accomplished through new 
measurements [87]. The issue of how the value of \V U &\ determined from neutron beta decay parameters r n 
and A compares with those from the superallowed decays and CKM unitarity can not be fully addressed 
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without a resolution of the value of t„. 

2.2.2 Prospects in Neutron Lifetime Measurements 

Several new experimental efforts to measure t„ are underway that include both refinement of existing ex- 
periments and techniques as well as development of new approaches. Most of the proposed experiments 
use confined UCN since the systematic effects are thought to be less challenging. The beam experiments 
are important because their set of systematic effects is almost completely orthogonal to those of the UCN 
confinement experiments. Agreement using the two techniques is needed for reliability of the value of the 
lifetime. 

One experiment proposes to use a low temperature Fomblin oil to reduce the UCN collisional losses . 
Recent data demonstrate a UCN reflection loss coefficient of 5 x 10 -6 per bounce when the vessel temperature 
is in the range of 105 K to 150 K. They intend to use this surface coating in an "accordion-like" storage vessel, 
thus allowing one to vary the trap volume while keeping the surface area and characteristics constant [88]. 
The collaboration expects to achieve a precision of 1 s. Another experiment will store UCN magnetically 
in vacuum using an arrangement of permanent magnets and superconducting solenoids and extract the 
protons electrostatically. The lifetime is measured by real-time detection of the decay protons and counting 
the integral number of neutrons using different storage times [89,90]. The collaboration anticipates that a 
measurement with 0.1s uncertainty is possible. There are also ideas to use a quadrupole trap [91] and to 
extract UCN into an Ioffc-type magnetic trap [92]. 

A third approach to measuring r„ exists that, in principle, avoids this set of systematic problems. The 
most natural way to measure exponential decay is to acquire an ensemble of radioactive species and register 
the decay products. One can then simply fit the resulting time spectrum for the slope of the exponential 
function. With advances in the density of stored ultracold neutrons and magnetic field technology, a number 
of groups propose to eliminate the losses from wall scattering in material bottles and trap the neutrons 
in a combination of magnetic field gradients and gravity [93]. In this way the neutrons never encounter 
a dissipative interaction during the measurement and the loss of neutrons from the trap should be caused 
only by neutron decay. This class of experiments that continuously monitor neutron decays must become 
statistically competitive to make a significant statement, but the prospects are encouraging. 

One experiment that has made significant progress uses ultracold neutrons that are magnetically confined 
in superfiuid 4 Hc [94]. The UCNs fill a magnetic trap through the inelastic scattering of 0.89 nm neutrons 
in superfiuid 4 He (the superthermal process). The decay electrons are registered via scintillations in the 
helium thus allowing one to directly fit for the exponential decay of the trapped neutrons. The first result 
is in agreement with the currently accepted value of the free neutron lifetime, but the statistical uncertainty 
is large (60s) [95]. Upgrades to the apparatus are in progress to increase the number of trapped neutrons. 
The collaboration anticipates that a statistical precision of a few seconds will be possible in the near future. 

2.2.3 Radiative Neutron Decay 

Beta decay of the neutron into a proton, electron, and antineutrino is accompanied by the emission of an 
inner-bremsstrahlung (IB) photon. While IB has been measured in nuclear beta decay and electron capture 
decays, it has never been observed in free neutron beta decay until very recently. Gaponov and Khafi- 
zov calculated the photon energy spectrum and branching ratio within a quantum electrodynamics (QED) 
framework [96-98], while Bernard et al. have calculated those properties as well as the photon polarization 
using heavy baryon chiral perturbation theory (HB%PT) including explicit A degrees of freedom [60,61]. 
The QED calculation takes into account electron IB while the HBxPT calculation includes all terms to 
order 1 /M (where M is the nucleon mass) including photon emission from the effective weak vertex. These 
additional terms contribute at the level of less than 0.5 % and create only a slight change in the photon 
spectrum and branching ratio. Both the photon energy spectrum and the photon polarization observables 
are dominated by electron IB. 

The experimental challenge was to distinguish the very low rate of radiative decay events in the large 
photon background of a neutron beam. The branching ratio above 15 keV is only about 3 x 10~ 3 , which (cou- 
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pled with the long neutron lifetime) makes the rate of detectable photons quite small. The first experiment 
to search for this decay mode resulted in an upper limit of 6.9 x 1CT 3 (at the 90 % confidence level) for the 
branching ratio of photons between 35 keV and 100 keV [99]. A more recent effort was able to definitively 
detect the radiative photons and measure the branching ratio at the 10% level [52]. In the experiment, 
photons with energies between 15keV and 340 keV were detected by a scintillating crystal coupled to an 
avalanche photodiode and were distinguished from uncorrelated background photons by coincidence with 
both the decay electron and proton [100]. Correlated background from external bremsstrahlung generated 
in the electron detector was estimated to contribute less than 3% of the radiative decay event rate. The 
branching ratio was measured to be (3.13 ± 0.34) x 10~ 3 and is consistent with the theoretical predictions 
of 2.85 x 10~ 3 in the same energy region. Only a small fraction of the solid angle for photon detection was 
utilized in the experiment. A new detector is under construction that should permit a precision measurement 
of the photon energy spectrum and the branching ratio at the few percent level. A measurement below the 
0.5 % level could reveal direct emission from the weak vertex. Furthermore, a measurement of the photon 
circular polarization could reveal information about the Dirac structure of the weak current [60,98]. Access 
to the radiative photon opens the possibility of new areas of investigation for neutron beta decay. 

2.3 Angular Correlation Experiments 
2.3.1 Spin-electron Asymmetry A 

With the neutron lifetime and one of the correlation coefficients a, A, or B, one can determine values for 
gA and gy. Because it has the greatest sensitivity to A and is more accessible experimentally, the spin- 
electron asymmetry A has been measured more frequently and with greater precision. Four independent 
measurements used in the PDG evaluation are not in good agreement with each other, so the PDG uses a 
weighted average for the central value and increases the overall uncertainty by a scale factor of 2.3 [57]. We 
discuss the two more recent measurements, one using a time projection chamber and one using an electron 
spectrometer, and the prospects for future improvement. 

In the experiment of Schrcckenbach et al. [101], a beam of polarized cold neutrons was surrounded by 
a time projection chamber (TPC). Decay electrons passed through the drift chamber and were incident 
on plastic scintillators. The drift chamber recorded the ionization tracks in three dimensions while the 
scintillator gave the electron energy and start signal for the drift chamber. The TPC provided good event 
identification and reduced gamma ray backgrounds. The result was A = —0.1160 ± 0.0015 [102]. The 
contributions to the overall uncertainty were roughly split between statistical and systematic uncertainties 
with the largest systematic contribution coming from the background subtraction. 

The PERKEO II experiment also used a beam of cold polarized neutrons, but the decay electrons were 
extracted using a superconducting magnet in a split pair configuration. The field was transverse to the beam, 
so neutrons passed through the spectrometer but electrons were guided by the field to one of two scintillator 
detectors on each end. This arrangement had the advantage of achieving a 4ir acceptance of electrons. They 
form an asymmetry from the electron spectra in the two detectors as a function of the electron energy; 
the difference in those quantities for the two detectors is directly related to the electron asymmetry. Their 
run produced A = — 0.1 178 ± 0.0007, where the main contributions to the uncertainty were in the neutron 
polarimetry, background subtraction, and electron detector response [72]. 

The next version of PERKEO will use the new ballistic supermirror guide at the ILL with four times 
the fluence rate [22] . They intend to use a new configuration of crossed supermirror polarizers to make the 
neutron polarization more uniform in phase space [103, 104]. The beam polarization can also be measured 
with a completely different method using an opaque 3 Hc spin filter [105]. They anticipate reducing the main 
correction and uncertainty in the polarization analysis from 1.1 % to less than 0.25 % with an uncertainty of 
0.1 % in that value. 

There arc several other efforts underway to perform independent measurements of the electron asymme- 
try. The UCNA collaboration has made progress toward measuring A using a superconducting solenoidal 
spectrometer [106]. UCN are produced in a solid deuterium moderator at LANSCE and transported to the 
spectrometer using diamond-coated guides. In the spectrometer, one produces highly polarized (> 99.9%) 
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neutrons by passing them through a 6 T magnetic field and into an open ended cylinder which increases the 
dwell time of the polarized UCN in the decay spectrometer. The decay electrons will be transported along 
the field lines to detectors at each end of the spectrometer. The detectors will initially consist of multiwire 
proportional counters backed by plastic scintillator, but other detectors are under development. Extensive 
measurements of electron backscattering on electron detector materials have been performed [107, 108] The 
collaboration believes that a 0.2 % measurement is possible with three weeks of running. 

Two other groups propose measuring the A coefficient at the 10~ 3 level. The detector designs allow the 
possibility of measuring other decay correlation coefficients with the same apparatus. A group at PNPI is 
working on a magnetic spectrometer to be used with a highly collimated cold neutron beam [109] . The field 
guides decay particles to an electron detector at one end and a proton detector at the other. Their neutron 
polarimeter agrees with 3 He-based spin filter methods at the «2x 10~ 3 level [110]. 

The abBA collaboration proposes to use an electromagnetic spectrometer that guides both decay elec- 
trons and protons to detectors at each end of the spectrometer [111]. The detector would be able to measure 
a, A, B, and the Fierz interference term b all in the same apparatus. The detectors would be large-area 
segmented silicon detectors with thin entrance windows that allow the detection of both the proton and 
electron. The ability to detect coincidences greatly suppresses backgrounds and allows the measurement of 
residual backgrounds. The magnetic field guides the decay products to conjugate points on the segmented 
Si detectors and provides 47r detection of both electrons and protons and suppression of backgrounds by use 
of coincidences. The apparatus is being designed for use at a pulsed spallation source to exploit background 
reduction and perform neutron polarimetry. The neutrons can be polarized by transmission through po- 
larized 3 He, whose spin-dependent absorption cross section possesses an accurately-known neutron energy 
dependence that can be exploited for accurate neutron polarization measurement [105,112-117]. 

One strategy to eliminate the need for absolute neutron polarization measurements is to measure both 
PA and PB simultaneously in the same apparatus (P is the neutron polarization) and take the ratio. From 
this ratio A can be calculated directly. Many of the planned experiments propose to implement this method, 
which was first performed in the experiment of Mostovoi [118] with the result A = —1.2686 ± 0.0046 

2.3.2 Spin-antineutrino Asymmetry B 

The electron asymmetry and antineutrino asymmetry provide complementary information. §a is equal to 
— 1 in the SM Lagrangian at the quark level but is renormalized in hadrons by the strong interaction. Since 
qa is nearly —1, A is close to zero, and B is near unity, B is not particularly sensitive to A but is more 
sensitive to certain non-SM extensions, such as extended left-right symmetric models. 

Left-right symmetric models, which are motivated in part to restore parity conservation at high energy 
scales, add a new right-handed charged gauge-boson Wi with mass M 2 and four new parameters to be 
constrained by experiments: a mixing angle £, S — (M1/M2) 2 , r g — gn/gh the ratio of the right- and left- 
handed gauge coupling strengths gn and gL, and Rk = Kid/Kni wnere R and L designate the right and left 
sectors. In the manifest left-right symmetric model (MLRM), r g = rjf = 1, and in the SM (5 = 0. 

The best direct search for W2 sets a lower limit of 652GeV/c 2 on M2 [119], and the mass limit on a 
right-handed vector boson comes from muon decay and is 406GeV/c 2 [120]. In the MLRM where there are 
only two parameters (( and 5), constraints from other systems are better than the neutron constraints. For 
the extended left-right model, however, neutron-derived constraints are complementary to the other searches. 
Another area in which to search for right-handed currents is the decay of the neutron into a hydrogen atom 
and antineutrino since one of the hyperfine levels of hydrogen in this decay mode cannot be populated unless 
right-handed currents are present [51]. The small branching ratio has precluded a search so far. 

In the last three decades, there have been only two new measurements of the antineutrino asymmetry. 
Since the antineutrino cannot be conveniently detected, its momentum was deduced from electron-proton 
coincidence measurements. Electrons from the decay of polarized neutrons were detected by plastic scintil- 
lators, and protons were detected by an assembly of two microchannel plates. From the electron energy and 
proton time-of-flight, one can reconstruct the antineutrino momentum. The first measurement was carried 
out at PNPI and produced a result of B = 0.9894 ± 0.0083 [121]. A second run at the ILL used largely 
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the same apparatus and measured B = 0.9821 ± 0.0040 [122], where the largest reduction in the overall 
uncertainty came from improved statistics. 

A recent measurement of B was performed using the PERKEO II apparatus [123]. Typically, one detects 
electron-proton coincidences using one detector for each particle. The PERKEO II measurement uses two 
detectors, one in each hemisphere of the detector, that can detect both electrons and protons. Electrons are 
detected using plastic scintillator, while the protons are accelerated on a thin carbon foil placed in front of 
the scintillator. The resulting secondary electrons are guided onto the electron detectors. This technique 
reduces systematics and increases the sensitivity to B. The result, B = 0.967 ± 0.006(stat) ± 0.010(sys) can 
be further improved in PERKEO III. The PERKEO II apparatus also lends itself to the measurement of 
the proton asymmetry coefficient C, and the first measurement of C has recently been reported [124]. 

2.3.3 Electron-antineutrino Asymmetry a 

Although the electron-antineutrino asymmetry a has approximately the same sensitivity to A as A, it is only 
known to 4%. Since 1978 [125], there has been only one new measurement. Unlike A, its measurement 
does not require polarized neutrons. The experimental difficulty lies with the energy measurement of the 
recoil protons, whose spectral shape is slightly distorted for nonzero a. A precision measurement of a would 
produce an independent measurement of A, an improved test of CKM unitarity, and model-independent tests 
of new physics. The values of a, A, and B can be related to the strength of hypothetical right-handed weak 
forces and scalar and tensor forces [126, 127], and it was recently shown that a precise comparison of a and 
A can place stringent limits on possible conserved-vector-current (CVC) violation and second class currents 
in neutron decay [128]. 

The most recent determination of a comes from measurements of the integrated energy spectrum of recoil 
protons stored in an ion trap [129, 130]. A collimated beam of cold neutrons passed though a proton trap 
consisting of annular electrodes coaxial with a magnetic field whose strength varied from 0.6 T to 4.3 T over 
the length of the trap. Protons created inside the volume were trapped, and those created in a high field 
region were adiabatically focused onto a mirror in the low field region. The trap was periodically emptied 
and the protons counted as a function of the mirror potential. The result of a = —0.1054 ± 0.0055 is in good 
agreement with the previous measurement and of comparable precision. 

The precision of a measurements must be improved to the level of A experiments to constrain A. There 
are two major efforts underway to improve the precision of a, aSPECT [131] and aCORN [132, 133]. aCORN 
relies on the measurement of an asymmetry in the coincidence detection of electrons and recoil proton that 
is proportional to a [134]. The asymmetry is formed by carefully restricting the phase space for the decay 
in a magnetic spectrometer so that decay events with parallel and antiparallel electron and antineutrino 
momenta are separated in the coincidence timing spectrum, a is directly proportional to the relative number 
of events, and there is no need for precise spectroscopy of the low energy protons. The experiment will 
be built and tested at the Low Energy Neutron Source (LENS) [135, 136] and then run at NIST where a 
measurement of approximately 1 % accuracy is feasible. 

In the aSPECT experiment, one again measures a proton energy spectrum as a function of a potential, 
similar to the idea used for the proton trap experiment. One increases the statistical power by completely 
separating the source part and the spectroscopy part of the apparatus. A cold neutron beam will pass 
through a region of strong, homogeneous magnetic field transverse to the beam. The decay protons with 
initial momentum component along the field direction will be directed toward a detector. Near the detector 
is a region of weaker magnetic field and electrostatic retardation potentials, and only those protons with 
sufficient energy to overcome the barrier continue on to the detector. Registering the protons as a function of 
the retardation potential gives the recoil proton spectrum, which one fits to extract a. This type of integral 
spectrometer is similar in spirit to the latest electron spectrometers used in tritium beta decay searches for 
the electron neutrino mass. The collaboration believes that a statistical uncertainty of approximately 0.25 % 
is achievable [137] and has commenced testing the apparatus at the ILL. Other ideas for measuring a have 
recently been proposed [138]. 
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3 SEARCHES FOR NONSTANDARD T AND B VIOLATION 



The physical origins of the observed CP violation in nature, first seen in the neutral kaon system [139], 
remain obscure. CP violation implies time-reversal symmetry T violation (and vice versa) through the CPT 
theorem. Recent experiments have reported measuring CP violation in the K° — ► 2n amplitudes [140, 141] 
and in the decays of the neutral S-mesons [142,143]. The Standard Model can accommodate the possibility 
of CP violation through a complex phase 5km in the CKM quark mixing matrix. To date there is no 
firm evidence against the possibility that the observed CP- violation effects are due to this phase [144], but 
the question remains whether or not there are sources of CP- violation other than 5km- There is indirect 
evidence for this possibility from cosmology; it appears that 5km is not sufficient to generate the baryon 
asymmetry of the universe in the Big Bang model. One area to probe for the existence of new CP-violating 
interactions is systems involving first-generation quarks and leptons for which the contribution from 5km 
is typically suppressed. Examples of observables of this kind are electric dipole moments of the neutron, 
leptons, and atoms [145] and T-odd correlations in leptonic and semileptonic decays and reactions and 
neutron spin-dependent transmission asymmetries. 

3.1 EDM Theoretical Framework 

The search for the neutron electric dipole moment addresses issues which lie at the heart of modern cosmology 
and particle physics. The current limit on the permanent EDM of the neutron represents one of the most 
sensitive null measurements in all of physics and has eliminated many theories and extensions to the Standard 
Model (see Figure |4j. The reader is directed to comprehensive reviews of EDM experiments in Refs. [47, 
146-148]. 

The energy of a neutral spin- 1/2 particle with an electric dipole moment d n in an electric field E is 
En = d n a ■ E, where a is the Pauli spin matrix. This expression is odd under T and P. The current 
experimental bound on the neutron EDM is d n < 0.63 x 10~ 25 e • cm (90 % CL) [57, 149]. In the SM, there 
are two sources of CP violation. One source is the complex phase 8km in the CKM matrix. The other 
source is a possible term in the QCD Lagrangian itself, the so-called 0-term 

£qcd = £qcd,8=o + Sr^GuvG^ , (9) 

6Z7T 

which explicitly violates CP symmetry because of the appearance of the product of the gluonic field operator 
G and its dual G. Since G couples to quarks but does not induce flavor change, d n is much more sensitive to 
9 than it is to 8km- Thus, measurement of d n determines an important parameter of the SM. Calculations 
have shown that d n ~ O(10~ 16 6>) e • cm [150-152]. 

The observed limit on d n allows one to conclude that 9 < 10~( 9±1) [153]. Since the natural scale is 
9 ~ 0(1), the very small value for 9 (known as the strong CP problem) requires an explanation. One 
attempt augments the SM by a global U(l) symmetry (referred to as the Peccei-Quinn symmetry), whose 
spontaneous breakdown gives rise to Goldstone bosons called axions [154,155]. The #-term is then essentially 
eliminated by the vacuum expectation value of the axion. No axions have yet been observed despite extensive 
searches. 

Since CP violation through the phase in the CKM matrix involves flavor mixing of higher generation 
quarks, d n is very small in the SM; calculations predict it to be 10~ 32 e • cm to 10~ 31 e • cm [156, 157], several 
orders of magnitude beyond the reach of any experiment being considered at present. Models of new physics, 
including left-right symmetric models, non-minimal models in the Higgs sector, and supersymmetric (SUSY) 
models, allow for CP violating mechanisms not found in the SM, including terms that do not change flavor. 
Searches for electric dipole moments in the neutron, leptons, and atoms, which are particularly insensitive 
to flavor-changing parameters, can therefore strongly constrain such models. 
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Figure 4: Plots showing the history of neutron EDM limit and some of the ranges for different theoretical 
predictions. (Plot courtesy of P. Harris, University of Sussex.) 
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3.1.1 Baryon Asymmetry 

Antimatter appears to be rare in the universe, and there is an asymmetry between the number of baryons 
and antibaryons. Although the SM possesses a nonperturbative mechanism to violate B, no experiments 
have seen B violation, and it is natural to speculate on the origin of the baryon asymmetry of the universe. 
There are two outstanding facts: baryons make up only 5% of the total energy density of the universe and 
the ratio of baryons to photons is very small. The ratio ns/n 7 = (6.1 ± 0.3) x 10~ 10 is known independently 
both from BBN and fluctuations in the microwave background [158]. Sakharov first raised the possibility of 
calculating the baryon asymmetry from basic principles [159]. He identified three criteria that, if satisfied 
simultaneously, will lead to a baryon asymmetry from an initial B = state: baryon number violation, CP 
violation, and departure from thermal equilibrium. One way to explain the asymmetry assumes that all 
three of these conditions were met at some very early time in the universe and that this physics will remain 
inaccessible to us, with the B asymmetry effectively an initial condition. However, the scenario of inflation 
in the early universe - which generates a flat universe, solves various cosmological problems, and generates 
a spectrum of primordial density fluctuations consistent with observation - would dilute any such early B 
asymmetry to a negligible level. In this case the B asymmetry must be regenerated through subsequent 
post-inflationary processes, and there is hope that it is calculable from first principles [160]. 

Although the SM contains processes that satisfy the first two conditions and the Big Bang certainly 
satisfies the third, it fails by many orders of magnitude in its estimate of the size of the baryon asymmetry. 
Grand unified theory (GUT) baryogenesis at T <~ 10 29 K corresponding to a mass scale on the order of 
10 16 GeV is disfavored by inflation. Electroweak baryogenesis [161, 162], which relics on a nonperturbative 
B — L-violating mechanism present in the SM due to nonperturbative electroweak fields [163] combined 
with CP violation and a departure from equilibrium at the electroweak phase transition, remains a viable 
possibility. Leptogenesis [164, 165] combined with B — L conserving processes to get the B asymmetry and 
the Affleck-Dine mechanism are favored speculations at the moment [166]. 

It appears that some physics beyond the SM, including new sources of CP violation that may lead to 
a measurable value for d n , must exist if the observed baryon asymmetry is to be understood. The minimal 
supersymmetric extension of the SM (MSSM) [167] can possess small values of the CP-violating phases 
(consistent with constraints from d n ) that generate the baryon asymmetry. Within the broad framework of 
non-minimal SUSY models, including GUTs, there are numerous new sources of CP violation in complex 
Yukawa couplings and other Higgs parameters that may have observable effects on the neutron EDM [168— 
170]. The present limit on the neutron EDM already severely constrains many SUSY models. The question 
of the relationship between electric dipolc moments and the baryon asymmetry generated at the electroweak 
scale is just beginning to be addressed in a quantitative manner [171] 

3.2 Electric Dipole Moment Experiments 

EDM experiments monitor the possible electric field dependence of the precession frequency of polarized 
neutrons using the well-known Ramsey interferometric technique of separated oscillatory fields. Static electric 
Eq and magnetic Bq fields are applied to the polarized neutrons. The neutron spin state is then governed by 
the Hamiltonian H = —jl ■ B n ± d n ■ E . An RF magnetic field of frequency cj a is applied to tilt the neutron 
polarization normal to Eo and Bo and it starts to precess with a frequency ujr. After a free precession time 
T a second tilt pulse in phase with the first is applied and the neutron polarization direction is proportional 
to (wr — Lo a )T. The Larmor precession frequency of the neutron depends on the direction of the applied field 
Eq relative to jl. An EDM would appear as a change in ujr as the electric field is reversed. 

The two most stringent limits on d n come from Altarev et al. at PNPI [172] and Harris et al. at 
ILL [149]. Both experiments used stored UCN. The PNPI apparatus contained two UCN storage chambers 
with oppositely-directed electric fields. A nonzero EDM would cause frequency shifts of opposite sign in each 
of the chambers, and some sources of magnetic field noise are suppressed with simultaneous measurements 
with both fields. Nevertheless, slowing varying magnetic fields remained a significant source of systematic 
uncertainty. In the experiment of Harris et al., a polarized 199 Hg comagnetometer occupying approximately 
the same volume as the neutrons was introduced into the storage volume to continuously monitor the 
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magnetic field. The comagnctometcr was essential in eliminating stray magnetic fields as a major source of 
systematic uncertainty. The experimental accuracy was limited by neutron counting statistics. 

There are ambitious efforts underway to improve the current neutron EDM limit by one to two orders 
of magnitude. All of the experiments attempt to increase the number of UCN, the observation time, and 
the size of the applied electric field. As the measurements become more sensitive more subtle systematic 
effects, such as a recently-discussed frequency shift associated with certain UCN trajectories in a magnetic 
field gradient [173,174], must be considered. The CryoEDM collaboration intends to produce UCN through 
the superthermal process and transport them to a separate measurement chamber containing superfluid 
helium. Liquid helium should allow electric field values that arc several times larger than used in past exper- 
iments, and the cryogcnically pure environment should permit longer UCN storage times. The collaboration 
proposes to use a multichamber spectrometer for compensation of field fluctuations by means of SQUID 
magnetometers. 

The nEDM collaboration proposes to search for the neutron EDM using a double chamber storage cell 
to suppress magnetic field fluctuations and allow one to extract d n from the simultaneous measurement of 
chambers with opposite electric field values. The magnetic field would be inferred from a set of laser optically 
pumped Cs magnetometers placed outside the storage cells. The UCN would come from the solid deuterium 
UCN source under construction at PSI. 

Another EDM experiment under development [148, 175] also proposes to increase the UCN density using 
downscattering of UCN in superfluid 4 He and exploit the large electric fields achievable in helium. The exper- 
iment will use polarized 3 He atoms as the comagnetometer in a bath of superfluid helium at a temperature of 
approximately 300 mK. The comagnetometer should sample the same fields as the neutrons, and historically 
this has been a very important check on potential systematic errors from magnetic field fluctuations. The 
strong spin dependence of the 3 He neutron absorption cross section allows the relative orientation of the 
neutron and 3 He spins to be continuously monitored through the intensity of the scintillation light in the 
helium in addition to SQUID magnctometry techniques. The polarized 3 He slowly depolarizes in the storage 
chamber and must be periodically replenished: the low temperature dynamics of 3 He motion in superfluid 
4 He makes this possible, and unpolarized test measurements have been performed [176,177]. One obtains 
d n by measuring the difference in the neutron and 3 He precession frequencies for the different orientation of 
the electric field, which can be measured optically in superfluid helium using the Kerr effect [178]. 

There are also preparations underway to search for the neutron EDM using dynamical diffraction from 
noncentrosymmetric perfect crystals. In dynamical diffraction the incident neutron plane wave state \k > 
is split as it enters the crystal into two coherent branches \k + > and \k- > with slightly different momenta 
and energies. The probability density of these two states in the crystal is concentrated along and in between 
the lattice planes, respectively. In noncentrosymmetric crystals the position of the electric field maxima can 
be displaced with respect to the nuclei. Therefore one of the branches can experience interplanar electric 
fields(~ 10 9 V/cm) which are orders of magnitude larger than can be achieved through application of external 
fields in vacuum or in isotropic media [179, 180]. The presence of a neutron EDM would produce an extra 
relative phase shift between the two interfering branches. Such experiments must contend with potentially 
large systematic effects such as those from neutron spin-orbit scattering from the atoms [181]. A number 
of experiments which investigate neutron optical issues relevant for an eventual EDM measurement of this 
type have been performed recently [182-188]. 

3.3 T-violation in Neutron Beta Decay 

With its small SM values of time-reversal violating obscrvablcs, neutron beta decay also provides a excellent 
laboratory in which to search for T violation. Lcptoquark, left-right symmetric models, and exotic fermion 
models can all lead to violations of time-reversal symmetry at potentially measurable levels [189]. One 
possible T-odd correlation in polarized neutron decay is Da n ■ (p e x p p ), where p p is the momentum of the 
recoil proton. The D coefficient is sensitive only to T-odd interactions with vector and axial vector currents. 
In a theory with such currents, the coefficients of the correlations depend on the magnitude and phase of 
A = |A|e-**. 
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D has T-even contributions from phase shifts due to pure Coulomb and weak magnetism scattering. The 
Coulomb term vanishes in lowest order in V-A theory [55], but scalar and tensor interactions could contribute. 
Fierz interference coefficient measurements [190, 191] can be used to limit this possible contribution to 
|£)EM| ^.8 x 10~ 5 )^. Interference between Coulomb scattering amplitudes and the weak magnetism 

amplitudes produces a final state effect of order E e 2 /p e m n . This weak magnetism effect is predicted to be 
j^VFMj _ i i x iq-5 [192]. References [144,189] summarizes the current constraints on D from analyses of 
data on other T-odd obscrvables for the SM and extensions. 

The EDM violates both T and P symmetries, whereas a D coefficient violates T but conserves P. 
This makes the two classes of experiments sensitive to different SM extensions. Although constraints on 
T-violating, P-conserving interactions can be derived from EDM measurements, these constraints may be 
model dependent [193], and EDM and neutron decay searches for T violation are complementary. 

3.4 D- and /^-coefficient Measurements 

In the last decade, there have been two major experimental efforts, EMIT and Trine, to improve the limit on 
the D coefficient in neutron decay. Each requires an intense, longitudinally polarized beam of cold neutrons 
around which one places proton and electron detectors alternating in azimuth. Coincidence data are collected 
in electron-proton pairs as a function of the neutron spin state to search for the triple correlation. 

In the emiT experiment, the detector consisted of four electron detectors and four proton detectors 
arranged octagonally around the neutron beam [194]. The octagonal geometry maximized the experiment's 
sensitivity to D by balancing the sine dependence of the cross product a n ■ (p e x p p ) with the large angles 
between the proton and electron momenta that are favored by kinematics. The decay protons drifted in a field 
free region before being focused by a 30 kV to 37 kV potential into an array of PIN diode detectors. With its 
maximum recoil energy of 750 eV, most of the protons arrived approximately 1 /is after the electrons. Detector 
pairs were grouped in the analysis to reduce potential systematic effects from neutron transverse polarization. 
The result from the first run of EMIT yielded an improved limit of D = [-0.6±1.2(stat)±0.5(sys)] x 10~ 3 [194]. 

Currently, the best constraint on D comes from the Trine collaboration, which reports D = [—2.8 ± 
6.4(stat)±3.0(sys)] x 10~ 4 [195]. They used two proton detectors and two electron detectors in a rectangular 
geometry. The proton detectors were comprised of arrays of thin- window, low- noise PIN diodes. The 
detectors were held at ground while the neutron beam was set to a potential of 25 kV by surrounding it 
with a high voltage electrode; the field was shaped to focus the decay protons onto the PIN arrays. The 
electrons were detected by plastic scintillators in coincidence with multi-wire proportional chambers. This 
coincidence provides reduction in the gamma-ray background rates and positional information on the decay, 
thus minimizing some sources of systematic uncertainty. 

The current PDG evaluation obtains a new value for the neutron D coefficient of (—4 ±6) x 10~ 4 , which 
constrains the phase of gA/gv to 180.06° ± 0.07° [57]. Neither experiment produced a statistically limited 
result, and both collaborations upgraded their detectors and performed second runs [196, 197]. In the near 
future it is reasonable to anticipate new results that will put a limit on D very near 10~ 4 . Although there 
have been discussions and ideas for experiments using UCNs, there are currently no concrete proposals to 
further improve the limit on D. 

Another T-odd correlation that may be present in neutron decay is the R correlation, R<j* n ■ (a e x p e ), 
where a n is the neutron spin and a e is the spin of the decay electron. A nonzero R requires the presence 
of scalar or tensor couplings and is sensitive to different SM extensions than D. An effort is underway at 
PSI to measure R in the neutron by measuring the neutron polarization and the momentum and transverse 
polarization of the decay electron at the level of 5 x 10~ 3 [198]. 

Neutron decay is a mixed Fermi and Gamov- Teller decay, so a measurement of R would produce a 
limit on both scalar and tensor T-odd couplings. The limit on R achieved in 8 Li Gamov-teller decay of 
R = (0.9 ± 2.2) x 10~ 3 now sets the most stringent limits for time-reversal violating tensor couplings in 
semileptonic weak decays, -0.022 < Im(C T + C' T )/C A < 0.017 [199]. 
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3.5 T Violation in Neutron Reactions 



T violation can lead to terms in the forward scattering amplitude for polarized neutrons in polarized or 
aligned targets of the form s ■ (k x fj and s ■ (k x T)(k ■ T) [200], where s is the neutron spin and / is 
the nuclear polarization. Since the enhancement mechanisms for parity violation in compound resonances of 
heavy nuclei are also applicable to T-odd interactions [201-203] (see Section^), it is possible for such searches 
to reach interesting levels of sensitivity. Although in principle T-odd observables in forward scattering are 
motion-reversal invariant and therefore not subject to final state effects, in practice the large spin dependence 
of the neutron- nucleus strong interaction in a polarized target can induce large potential sources of systematic 
errors which require careful study [204]. 

These systematic effects are smaller in aligned targets, and a search for the T-odd fivefold correlation 
in MeV polarized neutron transmission in an aligned holmium target has set the best direct limit on such 
interactions [205]. This P-even, T-odd correlation is especially interesting to search for, since there exists 
no renormalizable gauge theories with P-even T-odd tree-level gauge boson couplings between quarks [206] . 
Although EDM limits can also be used to constrain P-even T-odd interactions in many models, in general 
only direct measurements can set model- independent bounds [207] . A JINR-ITEP collaboration is preparing 
to perform a search for the P-even T-odd fivefold correlation with low energy neutrons on p-wave resonances 
using microwave-induced dynamical nuclear alignment to order the nuclei [208-210]. A shift in the zero 
crossing of the forward-backward asymmetry in unpolarized neutron capture on p-wave resonances has also 
been considered as a signal for P-even, T-odd interactions [211] but it is not a null test and although the 
statistical precision of the bound from 113 Cd is at the 10~ 4 level, one must accumulate a number of results 
either in different nuclei or in different resonances in the same nucleus in order to construct a statistically 
meaningful bound from such measurements [212]. 

Searches for the 3-fold P-odd T-odd correlation require a polarized target. Nuclei have been identified 
( 139 La, 131 Xe) that are polarizable in macroscopic quantities and possess large parity-odd asymmetries at 
low energy p-wave resonances [213-215]. The first steps toward such an experiment are in progress at 
KEK [216-219]. Polarized neutron capture on unpolarized targets can also set limits on P-odd, T-odd 
interactions [220], but because there are inelastic reactions, they are susceptible to initial and final state 
effects that can mimic T violation [221,222]. 

3.6 Neutron-antineutron Oscillations 

An observation of neutron-antineutron oscillations would constitute a discovery of fundamental impor- 
tance [223] . It requires a change of baryon number by 2 units and no change in lepton number and therefore 
must be mediated by an interaction outside the SM of particle physics. Among neutral mesons with suf- 
ficiently long lifetimes (kaons, B mesons) ) or neutral systems with flavor changing processes (neutrinos), 
oscillations are no longer a surprising phenomenon. The observation of oscillations in these systems has 
yielded information on aspects of physics (lepton number violation, T violation, neutrino mass) that are not 
accessible using less sensitive techniques. It is reasonable to hope that a search for oscillations in the neutron, 
the only neutral baryon which is sufficiently long-lived to conduct a practical experiment, may uncover new 
processes in nature. 

In the SM there are no (perturbative) renormalizable interactions one can write down which violate P, 
and any nonrenormalizable operator that can induce B violation must be naively suppressed by some heavy 
mass scale. The effective operator for neutron-antineutron oscillations involves a dimension 9 operator to 
change the 3 quarks in the neutron into 3 antiquarks and is suppressed by some mass scale to the 5th power. 
Some SM extensions lead to B violation by 2 units and not 1 unit. Examples include left-right symmetric 
models [224] with a local B-L symmetry needed to generate small Majorana neutrino masses by the seesaw 
mechanism, SUSY models with spontaneously broken B — L symmetry [225], and theories with compactified 
extra dimensions which attempt to solve hierarchy problems by introducing a much lower scale (TeV) for the 
onset of quantum gravity [226]. In these cases proton decay is unobservably small but neutron-antineutron 
oscillations can occur close to the present limit. A general analysis of all operators with scalar bilinears that 
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couple to two SM fermion fields uncovers operators which can only lead to ncutron-antineutron oscillations 
and not to proton decay [227]. 

The last experiment in the free neutron system at the ILL set an upper limit of 8.6 x 10 7 s (90 % confidence 
level) on the oscillation time [228]. Neutrons from the ILL cold source were extracted using a diverging 
neutron guide that compressed the transverse phase space of the neutron beam. The neutrons traversed 
a magnetically-shielded vacuum chamber of about 100 meters before striking a carbon foil surrounded by 
an antincutron detector. The residual magnetic field was measured in-situ by polarizing the neutrons and 
looking for the rotation of the plane of polarization using spin-echo techniques. No antineutron events were 
observed. Translated into a mass scale, this limit excludes mass scales for the effective operator that induces 
oscillations below ~ 100 TeV. A similar indirect limit is set by the absence of evidence for spontaneous 
neutron-antincutron oscillations in nuclei in large underground detectors built for proton decay and neutrino 
oscillation studies [229]. Although there has been some discussion of possible strategies to improve on the 
bounds from direct searches using cold and ultracold neutrons [230,231], there are no new free neutron- 
antincutron oscillation searches underway. 

The neutron and antineutron are measured to have the same mass to a precision of only 10~ 4 . If CPT 
symmetry is somehow violated and the neutron and antineutron mass are different, then neutron-antineutron 
oscillations are greatly suppressed [232]. 

There are many other modes one might imagine for £?-violating processes involving neutrons. Of the 
experimental limits for the many possible processes, the poorest limits are those for so-called "invisible" 
modes of neutron decay, favored in some models [233] in which one or two neutrons decay to neutrinos. 
If such processes were to occur in nuclei, they would leave neutron holes whose deexcitation would leave 
characteristic signatures. Recent data from underground detectors built to observe neutrino oscillations 
place much stronger constraints on such n and nn decay modes [234-236]. 

4 NEUTRON-NUCLEON WEAK INTERACTIONS 

The most obvious consequence of the weak interaction for neutrons is that it makes neutrons unstable. In 
addition to the coupling of quarks to leptons that allows neutron to decay, electroweak theory also predicts 
(and experiments confirm) that there are weak interactions between the quarks in the neutron with couplings 
comparable in size to those involved in neutron decay. The weak nucleon-nucleon (NN) interaction is a unique 
probe of strongly interacting systems. This section presents an overview of the importance of NN interactions 
for QCD and the status of the experimental efforts. Reviews of various aspects of the field can be found in 
References [202,237-241]. 

4.1 Overview 

The dynamics of the quarks in the nucleon are dominated by momentum transfers that are less than that set 
by the QCD scale of 1 GeV/c. In this regime QCD becomes so strong that quarks are permanently confined, 
and therefore the quark-quark weak interactions appear through the NN weak interactions that they induce. 
At these energies quark-quark weak amplitudes are of order 10~ 7 of strong amplitudes primarily because of 
the short range of the quark weak interactions through W and Z exchange. 

Assuming that it is correctly described by the electroweak theory at low energy scales, the quark-quark 
weak interaction can be used as an internal probe of strongly interacting systems. Collider measurements 
have verified the SM predictions for quark-quark weak couplings for large momentum transfers at the 10 % 
level [242]. The weak amplitudes are too weak to significantly affect the strong dynamics or to excite the 
system, and therefore it probes quark-quark correlations in the QCD ground state. The effects of the quark- 
quark weak interaction can be isolated from the strong interaction using parity violation. The short range of 
the quark-quark weak interaction and its ability to violate parity make it visible and sensitive to interesting 
aspects of strongly interacting systems, as seen in four cases. 

1) The ground state of the strongly interacting limit of QCD is a problem of fundamental importance. 
Although the dynamics that lead to the spontaneous breakdown of chiral symmetry in QCD are not yet 
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understood, one of the leading models assumes the importance of fluctuating nonperturbative gluon field 
configurations called instantons [243] . They induce four-quark vertices that flip the quark helicity and localize 
the quark wave function through a mechanism similar to Anderson localization of electrons in disordered 
metals [244,245]. Some aspects of QCD spectroscopy and the high density limit can be understood by 
assuming that quark-diquark configurations in the nucleon are important [246]. The mechanism for the 
phenomenon of color superconductivity in the high density limit of QCD consists of a BCS-like condensation 
of diquarks [247] . The quark-quark weak interaction in the nucleon in the low energy limit induces four-quark 
operators with a known spin and flavor dependence whose relative sizes are in principle sensitive to these 
and other correlation phenomena in the ground state of QCD. 

2) With experimental information on the low energy parity- violating (PV) partial waves in the NN system, 
there is a chance to understand quantitatively for the first time the extensive observations performed in many 
systems of parity-violating phenomena in nuclei [248] . Nuclear parity violation is linearly sensitive to small 
components of the nuclear wavefunction since successive shell model levels alternate in parity, and parity- 
odd operators directly connect adjacent shells [249]. Ideas from quantum chaos [250] and nuclear statistical 
spectroscopy have been used to analyze parity violation in neutron reactions in heavy nuclei in terms of 
the effective isovector and isoscalar weak NN interaction, and knowledge of PV in the NN system would 
allow a quantitative test of the predictive power of these ideas [251,252]. The matrix elements for weak NN 
interactions in nuclei also bear many similarities to the types of matrix elements that must be calculated to 
interpret limits on neutrino masses from double beta decay searches [253]. 

3) In atoms, the effect of NN parity violation was seen for the first time in 133 Cs [254] through its 
contribution to the anapolc moment of the nucleus, which is an axial vector coupling of the photon to the 
nucleus induced in heavy nuclei mainly by the PV NN interaction [255,256]. Anapolc moment measurements 
in other atoms are possible, and experiments are in progress [257]. In the heavy nuclei for which the anapolc 
moment is a well-defined observable, the main contribution comes from PV admixtures in the nuclear ground 
state wave function [258]. In electron scattering from nucleons, PV effects are sensitive to both Z exchange 
between the electron and the quarks in the nucleon as well as the coupling of the virtual photon to the axial 
current from PV interactions among the quarks in the nucleon. As the precision of PV electron scattering 
measurements improves in an effort to determine the strange magnetic moment of the nucleon [259-262], 
it may be important to know enough about the weak NN interaction to reliably extract the information of 
interest. 

4) The NN weak interaction is also the only practical way to study quark-quark neutral currents at low 
energy. The neutral weak current conserves quark flavor to high accuracy in the standard electroweak model 
(due to the GIM mechanism) . It is not seen at all in the well-studied strangeness-changing nonleptonic weak 
decays. We know nothing experimentally about how QCD modifies weak neutral currents. 

There are theoretical difficulties in trying to relate the underlying electroweak currents to low-energy 
observables in the strongly interacting regime of QCD. One expects the strong repulsion in the NN interaction 
to keep the nucleons sufficiently separated for a direct exchange of W and Z bosons between quarks to 
represent an accurate dynamical mechanism. If one knew weak NN couplings from experiment, they could 
be used to interpret parity violation effects in nuclei. The current approach is to split the problem into two 
parts. The first step maps QCD to an effective theory expressed in terms of the important degrees of freedom 
of low energy QCD, mesons and nucleons. In this process, the effects of quark-quark weak currents appear 
as parity- violating meson- nucleon couplings [263]. The second step uses this effective theory to calculate 
electroweak effects in the NN interaction in terms of weak couplings. The couplings themselves also become 
challenging targets for calculation from the Standard Model. 

4.2 Theoretical Description 

In the work of Desplanques, Donoghue, and Holstein (DDH) [263], the authors used a valence quark model in 
combination with SU(6) symmetry relations and data on hyperon decays to produce a range of predictions for 
effective PV meson-nucleon couplings from the SM. At low energy, the weak interaction between nucleons in 
this approach is parameterized by the weak pion coupling constant f n , and six other meson coupling denoted 
as hp, h p , hp 1 , h 2 p , h^, and h^, where the subscript denotes the exchange meson and the superscript indicates 
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the isospin change. Due to uncertainties in the effects of strong QCD, the range of predictions is rather 
broad. For the weak pion coupling, neutral currents should play a dominant role. Another strategy is to 
perform a systematic analysis of the weak NN interaction using an EFT approach to classify the interaction 
in a manner that is consistent with the symmetries of QCD and does not assume any specific dynamical 
mechanism. Such an EFT approach has recently appeared [264]. Preparations have also been made for an 
eventual calculation of the weak NN interaction vertices using lattice gauge theory in the partially quenched 
approximation [265]. 

At the low energies accessible with cold neutrons with k n R s t r ong << 1, parity-odd effects in the two- 
nucleon system can be parameterized by the five independent amplitudes for S — P transitions involving the 
following nucleons and isospin exchanges: 1 So — ► 3 Pq (p-p, p-n, n-n, AI = 0, 1, 2), 3 Si — > 1 P\ (n-p, AI = 0), 
and 3 Si — > 3 Pi (n-p, AI = 1). Thus, from the point of view of a phenomenological description of the weak 
NN interaction, at least five independent experiments are required. The PV longitudinal analyzing power in 
p-p scattering, which determines a linear combination of the 1 Sq — ► 3 pj amplitudes, has been measured at 
15MeV and 45MeV in several experiments with consistent results [266-269] and remains the only nonzero 
observation of parity violation in the pure NN system dominated by p-waves. 

Parity violation in few-nucleon systems should be cleanly interpretable in terms of the NN weak in- 
teraction due to recent theoretical and computational advances for the strong interaction in few nucleon 
systems [270]. Weak effects can be included as a perturbation. These calculations have recently been done 
for n-p and p-p parity violation [271-273] and can be done in principle for all light nuclei [274,275]. It is 
also possible that these microscopic calculations can be applied to systems with somewhat larger A, such as 
10 B and 6 Li, where measurements of P-odd observables with low energy neutrons have reached interesting 
levels of sensitivity [276,277]. 

The longest-range part of the interaction is dominated by the weak pion- nucleon coupling constant f n . f n 
has been calculated using QCD sum rules [278,279] and in a SU(3) Skyrme model [280], and a calculation in 
a nonlinear chiral quark model is in progress [281]. Measurements of the circular polarization of photons in 
the decay of 18 F [282,283] provide a value for / w that is considerably smaller than the DDH best value though 
still within the reasonable range. A precision atomic physics measurement of the 133 Cs hyperfine structure 
(anapole moment) has been analyzed to give a constraint on f n and the combination (h° p +0.6h^). This 
result would seem to favor a value for f n that is larger than the 18 F result. Figure [5] presents an exclusion 
plot that summarizes the current situation. Within the DDH one-boson exchange model of the weak NN 
interaction, all parity-odd amplitudes are actually products of strong and weak couplings, and changing the 
strong couplings to those from the Bonn potential can relax some of the observed discrepancies [284] 

4.3 Parity-Odd Neutron Spin Rotation and Capture Gamma Asymmetries 

There are a few general statements that apply to the low energy weak interactions of neutrons with low A 
nuclei. In the absence of resonances, it can be shown that the PV helicity dependence of the total cross 
section vanishes if only elastic scattering is present and that both the PV neutron spin rotation and the PV 
helicity dependence of the total cross section with inelastic channels are constant in the limit of zero neutron 
energy [200]. These results depend only on the requirement for parity violation in an S — ► P transition 
amplitude involving two-body channels. The two practical classes of neutron experiments are PV neutron 
spin rotation and PV gamma asymmetries. 

The NPDGamma experiment will measure the parity-violating directional gamma ray asymmetry A 1 in 
the capture of polarized neutrons on protons [286,287]. The unique feature of this observable is that it is 
sensitive to the weak pion coupling f„, A 1 — — 0.11/„- [288-290]. The recently-commissioned beamline at 
LANSCE delivers pulsed cold neutrons to the apparatus, where they are polarized by transmission through 
a large-volume polarized 3 He spin filter and are transported to a liquid parahydrogen target. A resonant 
RF spin flipper reverses the direction of the neutron spin on successive beam pulses using a sequence that 
minimizes susceptibility to some systematics. The 2.2 MeV gamma rays from the capture reaction are 
detected in an array of CsI(Tl) scintillators read out by vacuum photodiodes operated in current mode and 
coupled via low-noise I-V preamplifiers to transient digitizers. The current-mode Csl array possesses an 
intrinsic noise two orders of magnitude smaller than the shot noise from the gamma signal and has been 
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Figure 5: Constraints on linear combinations of isoscalar and isovector nucleon-nucleon weak meson cou- 
plings [285]. The bottom graphs show constraints from measurements of anapole moments in 133 Cs and 
205 T1 and from measurements in p-p and p- 4 He scattering and from 18 F and 19 F gamma decay. The top 
graph shows the anticipated constraints from proposed measurements of the P-odd asymmetry p(n,7)d to 
5 x 1CP 9 accuracy and the P-odd neutron spin rotation in 4 He to 2 x 10~ 7 rad/m accuracy. In each plot, 
the box indicates the DDH reasonable range for the couplings. 



shown in offline tests to possess no false instrumental asymmetries at the 5 x 10~ 9 level [291]. The pulsed 
beam enables the neutron energy to be determined by time-of- flight, which is an important advantage for 
diagnosing and reducing many types of systematic uncertainty. This apparatus has been used to conduct 
measurements of parity violation in several medium and heavy nuclei [292]. 

Another experiment in preparation is a search for parity violation in neutron spin rotation in liquid 
4 Hc [293]. A transverse rotation of the neutron spin vector about its momentum manifestly violates par- 
ity [294] and can be viewed from a neutron optical point of view as due to a helicity-dependent neutron 
index of refraction. For 4 He, the calculated PV neutron spin rotation in terms of weak couplings is [295] 

4> = {0.97 U + 0.32h° p - 0.11/iJ + 0.22/i° - 0.22/£) x 10" 6 rad/m (10) 

To measure the small parity-odd rotation, a neutron polarizer-analyzer pair is used to transmit only the 
component of a polarized beam that rotates as it traverses the target. The challenge is to distinguish small 
PV rotations from rotations that arise from residual magnetic fields. The first measurement achieved a 
sensitivity of 14 x 10~ 7 rad/m at NIST [296], and no systematic effects were seen at the 2 x 10~ 7 rad/m level. 
There are four plausible experiments which employ beams of cold neutrons and involve targets with 
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A < 5: measurement of the PNC gamma asymmetries in p(n,7)d and in d(n,7)t and of the PNC neutron 
spin rotations in 4 He and H. Successful measurements in all of these systems, in combination with existing 
measurements, would have a strong impact on the knowledge of the NN weak interaction [297]. 

4.4 Test of Statistical Theories for Heavy Nuclei Matrix Elements 

One might assume that a quantitative treatment of NN parity violation in neutron reactions with heavy nuclei 
would not be feasible. A low energy compound nuclear resonance expressed in terms of a Fock space basis in 
a shell model might possess a million components with essentially unknown coefficients, and the calculation 
of a parity-odd effect would involve a matrix element between such a state and another equally complicated 
state. However, one can imagine a theoretical approach which exploits the large number of essentially 
unknown coefficients in such complicated states. If we assume that it is possible to treat these components 
as random variables, one can devise statistical techniques to calculate, not the value of a particular parity- 
odd observable, but the width of the distribution of expected values. A similar strategy has been used to 
understand other properties of complicated compound nuclear states. The distribution of energy spacings 
has been known for a long time to obey a Porter-Thomas distribution [298] in agreement with the predictions 
of random matrix theory, and statistical approaches have been used to understand isospin violation in heavy 
nuclei [299] . Statistical analyses have been applied to an extensive series of measurements of parity violation 
in heavy nuclei performed mainly at Dubna, KEK, and LANSCE [300,301]. 

The TRIPLE collaboration at LANSCE measured 75 statistically significant PV asymmetries in several 
compound nuclear resonances in heavy nuclei. In the case of parity violation in compound resonances in 
neutron-nucleus reactions there are amplification mechanisms which can enhance parity-odd observables by 
factors as large as 10 5 . These amplification mechanisms, which are interesting phenomena in themselves, 
depend in an essential way on the complexity of the states involved and the reaction dynamics. Part of the 
amplification comes from the decrease in the spacing between levels as the number of nuclcons increases, 
which brings opposite parity states closer together and increases their weak mixing amplitudes [302,303], and 
for low energy neutron reactions in heavy nuclei it leads to a generic amplification of order 10 2 in parity-odd 
amplitudes. In addition, for low energy neutron-nucleus interactions the resonances are mainly I = and 
I = 1 with the scattering amplitudes in s-wave resonances larger than for p-wave resonances by a factor 
of order 10 2 to 10 3 . At an energy close to a p-wave resonance, the weak interaction mixes in an s-wave 
component that is typically much larger, and this factor also amplifies the asymmetry. These amplification 
mechanisms were predicted theoretically [304] before they were measured [305,306]. 

A basic tenet of the statistical approach is that there should be, on average, equal numbers of negative 
and positive PV asymmetries in a given nucleus. This appears to be satisfied in all nuclei measured so far 
except for 232 Th; below E = 250 eV all ten of the parity-odd asymmetries in this nucleus have the same 
sign. At present this result is ascribed to some poorly understood nuclear structure effect specific to 232 Th. 
Its observation illustrates the potential for the use of NN parity violation to discover unsuspected nuclear 
structure effects in complicated many body systems. In another recent example, a series of experiments 
that first observed parity-odd neutron spin rotation in natural Pb [307] and subsequent experiments on 
spin rotation in separated lead isotopes [308] led to a prediction of the first-known subthreshold p-wave 
resonance in neutron- nucleus interactions in the isotope with the large parity violation [309]. Although this 
isotope was thought to be 204 Pb according to the results of a spin rotation experiment at HMI [310], a 
subsequent neutron radiative capture experiment seems to indicate that 207 Pb possesses the subthreshold 
p-wave resonance [311]. 

The values of the weak matrix elements determined by the statistical analysis varied in the range 0.5 meV 
to 3.0 meV, in rough agreement with theory. The accuracy of this analysis was improved through new mea- 
surements of the required spectroscopic information on compound nuclear levels [312], and more spectroscopic 
information on the appropriate nuclei would be helpful. Theoretical calculations to use this data to extract 
the weak isoscalar and isovector couplings in 238 U obtain results in qualitative agreement with DDH expec- 
tations [252] . If the weak NN couplings were known, we should be able to use the results to see if there is 
any evidence for nuclear medium effects, which in principle should exist. Measurements of the mean square 
weak matrix elements in lower A regions would obviously be helpful. 
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4.4.1 Parity-odd and Time Reversal-odd Correlations in Neutron-induced Ternary Fission 

Another example where symmetry violation in neutron-induced reactions has led to progress in the under- 
standing of many-body nuclear dynamics is fission. P-odd effects in binary fission induced by polarized 
neutron capture have been observed for a long time [313,314]. Although one might expect that a treatment 
of parity violation in nuclear fission would be even more difficult than for compound nuclear resonances in 
heavy nuclei, there is a compelling understanding of parity-odd asymmetries observed in fission after capture 
by polarized neutrons [313-315] based on interference of amplitudes of opposite parity from parity doublets 
in the relatively small number of open tunneling channels through the cold pear-shaped transition state. 
Since this interference occurs among these small number of fission channels in the initial state near the 
saddle point, it can survive the inevitable averaging over the enormous number of final states later produced 
in the rupture. 

In the case of ternary fission, where a third light charged particle (usually an alpha) is emitted in addition 
to the two main fragments, recent parity violation measurements have given support to a specific mechanism 
for the emission of the ternary particle [316]. Consider two generic mechanisms: the simultaneous emission 
of the three particles (three-body compound nucleus decay) and "double neck rupture" in which the ternary 
particle is emitted after the first rupture of the neck from its remnants. In the first case, because all three 
objects originate from the same system where the dominant parity violation comes from the mixing of 
opposite parity compound nuclear states, one expects all of the PV asymmetries in various channels to be 
about the same size. In the second case, however, the mechanism for the emission of the ternary particle does 
not possess the same intermediate states that are known to exist in binary fission, and upon averaging over 
the large number of fragment states one would expect the parity-odd correlations that involve the ternary 
particle to be much smaller. This is what is observed experimentally in 233 U [317-321]. Furthermore, the 
parity-odd asymmetries of the two large fragments were seen to be independent of the energy of the ternary 
particle. Since different ternary particle energies are presumably coming from different Coulomb repulsion 
effects from different shapes of the neck, this independence would also seem to indicate that the parity-odd 
asymmetry is established before the scission process. 

In ternary fission one can also look for a formally T-odd triple correlation between the momenta of the light 
fragment and ternary particle and the neutron spin. This measurement has recently been done [322,323] and 
a large nonzero effect of order 10 -3 was seen in both 233 U and 235 U for both alphas and tritons as the ternary 
products. It is believed that in the ternary fission system this correlation is due to a final state effect and 
not to a fundamental source of T violation. The fact that the size of the observed triple correlation depends 
on the ternary particle energy also suggests that a final state effect is responsible. One model [324-326] 
can reproduce the order of magnitude of the effect if one assumes that the projection of the orbital angular 
momentum of the recoiling ternary particle changes the spin projections and therefore the level densities 
of the larger fragments. If the emission probabilities of the ternary particle are proportional to these level 
densities and the angular momentum of the initial system is correlated with the neutron polarization, this 
mechanism can generate a nonzero triple correlation. Semiclassically this can be viewed in terms of the 
Coriolis interaction of the ternary particle with the rotating compound nucleus [327]. Other recent work 
asserts that the large T-odd correlation in ternary fission is evidence for the simultaneous production of the 
three particles [328]. Future work will attempt to confirm this mechanism in plutonium. 

5 LOW ENERY QCD TESTS 

One of the long-term goals of strong interaction physics is to see how the properties of nucleons and nuclei 
follow from QCD. For nuclei the first step in this process is to see if one can start from QCD and calculate the 
well-measured NN strong interaction scattering amplitudes and the properties of the deuteron. Over the last 
decade a number of theoretical developments have started to show the outlines of how this connection between 
QCD and nuclear physics can be made. In this section, we discuss some of the theoretical developments 
in few nucleon systems along with the several precision scattering length experiments. The status of two 
fundamental properties of the neutron, its polarizability and the neutron-electron scattering length, are also 
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discussed. 



5.1 Theoretical Developments in Few Nucleon Systems and the Connection to 
QCD 

Based on a suggestion by Weinberg [329], one strategy to develop an effective field theory [330] for QCD 
that is valid in the low energy limit relevant for nuclei and incorporates the most important low energy 
symmetry of QCD is through chiral symmetry. It is recognized that this alone is not enough because some of 
the important energy scales of nuclear physics, such as the deuteron binding energy and its correspondingly 
large low energy scattering lengths, seem to be the result of a delicate cancelation between competing 
effects which will need more than chiral symmetry alone to understand. Recent efforts to understand the 
emergence of smaller energy scales in nuclear physics not set by chiral dynamics, such as the deuteron binding 
energy, have led to interesting suggestions that the low energy limit of QCD is not described by the usual 
renormalization group fixed point but rather is close to a limit cycle which can be reached by a small tuning 
of the values of the current quark masses in the Lagrangian [331,332]. 

Recently, significant insight into certain features of fcw-nuclcon systems has come from the EFT approach 
based on the chiral symmetry of QCD [333-335] . The value of the EFT approach is that it is a well-defined 
field theoretical procedure for the systematic construction of a low energy Lagrangian consistent with the 
symmetries of QCD. To a given level of accuracy the Lagrangian contains all possible terms consistent with 
symmetries with a number of arbitrary coefficients which, once fixed by experiment, can be used to calculate 
other observables. EFTs based on the chiral symmetry of QCD have been used to develop an understanding 
for the relative sizes of many quantities in nuclear physics, such as that of nuclear N-body forces [336] and in 
particular the nuclear 3-body force (3N), which is the subject of much activity. Although it is well understood 
that 3N forces must exist with a weaker strength and shorter range than the NN force, little else is known. 

EFT has been used to solve the two and three nucleon problems with short-range interactions in a 
systematic expansion of the small momentum region set by kb < 1, where k is the momentum transfer 
and b is the scattering length [265,335]. For the two-body system, EFT is equivalent to effective range 
theory and reproduces its well-known results for NN forces [337-339]. The chiral EFT expansion does not 
require the introduction of an operator corresponding to a 3N force until next-to-next-to leading order in 
the expansion, and at this order it requires only two low energy constants [340,341]. There have also been 
significant advances in other approaches to the computation of the properties of few-body nuclei with modern 
potentials [342] such as the AV18 potential [270,343], which includes electromagnetic terms and terms to 
account for charge-independence breaking and charge symmetry breaking. 

All of these developments show that precision measurements of low energy strong interaction properties, 
such as the zero energy scattering lengths and electromagnetic properties of small A nuclei, are becoming 
more important for strong interaction physics both as precise data that can be used to fix parameters in the 
EFT expansion and also as new targets for theoretical prediction. It is possible now to envision the accurate 
calculation of low energy neutron scattering lengths for systems with A > 2. 

5.2 Precision Scattering Length Measurements Using Interferometric Methods 

In parallel with these theoretical developments two interferometric methods have been perfected to allow 
high precision measurements of n-A scattering lengths. One is neutron interferomctry using diffraction from 
perfect silicon crystals [10], which measures the coherent scattering length. The other is pseudomagnctic 
precession of polarized neutrons in a polarized nuclear target, which measures the incoherent scattering 
length. Together these two measurements can be used to determine the scattering lengths in both channels. 

Neutron interferometry can be used to measure the phase shift caused by neutron propagation in the 
optical potential of a medium. For a target placed in one arm of an interferometer, the phase shift is given by 
the expression cf> — bNDX, where N is the target density, D is the thickness of the sample, A is the neutron 
de Broglie wavelength, and b is the coherent scattering length. Corrections to the relation between the phase 
shift and the coherent scattering length from multiple scattering effects in the theory of neutron optics are 
small and calculable [344]. High absolute accuracy in the determination of N,D, and A are required but 
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possible at the 10 4 level. Recent results from measurements at the NIST Neutron Interferometry and Optics 
Facility yielded the coherent scattering lengths b np = (-3.738 ± 0.002) fm, b nd = (6.665 ± 0.004) fm [345,346], 
and b n 3 He — (5.857 ± 0.007) fm [347] (see also [348]). They showed that almost all existing theoretical 
calculations of the n-d and n- 3 He coherent scattering lengths are in disagreement with experiment and that 
the accuracy of present measurements is sensitive to such effects as nuclear three-body forces and charge 
symmetry-breaking [349,350]. 

If the neutron- nucleus interaction is spin dependent, a polarized neutron moving through a polarized 
medium possesses a contribution to the forward scattering amplitude proportional to (b + — b-)a n ■ I where 
b + and &_ are the scattering lengths in the two channels, a n is the neutron spin, and I is the nuclear 
polarization. The angle of the polarization of a neutron initially polarized normal to the target polarization 
precesses as it moves through the medium. This phenomenon is referred to as nuclear pseudomagnetic 
precession [351] and has been used to measure scattering length differences in many nuclei. Recently, a high- 
precision measurement of this precession angle was performed in polarized 3 He using a neutron spin-echo 
spectrometer at the ILL [352]. A new experiment to determine the spin-dependence of the n-d scattering 
length is in preparation at PSI [353]. 

In combination with the n-d coherent scattering length determined by neutron interferometry, the PSI 
experiment should determine both n-d scattering lengths to 10~ 3 accuracy. The 2 Si/2 scattering length in 
the n-d system is especially interesting. The quartet s-wave scattering length ( 4 S 3 / 2 ) can be unambiguously 
determined from current theory. Because the three nuclcons in this channel exist in a spin-symmetric 
state, and hence have an antisymmetric space-isospin wavefunction, the scattering in this state is completely 
determined by the long range part of the triplet s-wave NN interaction in the n-p channel, i.e. by n-p 
scattering and the properties of the deuteron. By contrast the Pauli principle does not deter the doublet 
channel from exploring the shorter-range components of the NN interaction, where 3N forces should appear. 

Perhaps the single most interesting scattering length to measure is the neutron-neutron scattering length 
b nn . No direct measurements exist. An experiment to determine b nn by viewing a high-density neutron gas 
near the core of a pulsed reactor that produces an extremely high instantaneous neutron density and mea- 
suring a quadratic dependence of the neutron fluence on source power is currently being designed [354,355]. 
An experiment to let the neutrons in an extracted beam scatter from each other has been considered [356] . A 
precision (a few percent) measurement of b nn would be valuable, because it could give new information on the 
size of charge symmetry breaking effects of the strong interaction. Charge symmetry is the n-p interchange 
symmetry, a discrete subset of isospin symmetry corresponding to a 180 degree rotation in isospin space. A 
calculation of b nn using low energy effective field theories of QCD would seem to be out of reach at present 
since there are so few constraints on the nature of charge symmetry breaking in the strong interaction from 
other systems. EFT analyses to extract b nn from neutron-neutron final state effects in few body reactions 
are possible, and an EFT analysis of the n~d^ nn-y reaction to extract b nn has recently appeared [357]. 

It is interesting that both the n-d and n- 3 He scattering lengths are now of interest as input for few 
body calculations in other systems. The n-d coherent scattering length has been used recently to improve 
the precision of the determination of low energy constants in a chiral effective field theory for few nucleon 
systems [358]. The improvements in the n- 3 He scattering lengths are important input for a recent EFT 
calculation of the small cross section (but relevant in astrophysics) for the process n+ 3 He — > 4 He+7 [359] that 
was previously measured at the ILL [360]. These EFT calculations precisely reproduce radiative processes 
in few body systems such as the n +p — ► D + 7 total cross section [361]. A number of precise measurements 
have been performed on radiative processes in neutron capture on few body nuclei such as the total cross 
section for n + D — > 3 H + 7 and n + 3 He — > 4 He + 7 and the circular polarization inn + p^D + 7 [362] 
and n + D — > 3 H + 7 [363], and there is an upper bound on the gamma asymmetry in polarized neutron 
capture on polarized protons ri + p^ D + 7 [364]. As an example of the present level of sophistication of the 
theoretical calculations in this field one can see in where two and even three-body electromagnetic current 
operators are needed to reproduce the measurements of radiative capture in the n-D system [365]. EFT 
calculations, which are in principle traceable to QCD, can now be performed for these and other low energy 
observables, and it would be useful to perform more precise measurements in these systems to test the new 
predictive power of these approaches. 
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5.3 Neutron-electron Interaction 



The neutrality of the neutron is not a consequence of the minimal 3-generation Standard Model, and in 
principle it can possess a nonzero electric charge proportional to the difference in (conserved) lepton numbers 
of different generations [366] . The stringent upper limit on the electric charge of the neutron from laboratory 
experiments [367], q n < 10~ 21 e, strongly favors SM extensions, such as Majorana neutrino mass terms, 
which break family lepton number symmetries. One can imagine improving this direct limit using ultracold 
neutrons [368] although indirect limits from astrophysical arguments are much more stringent. There is also 
a stringent limit on a possible magnetic monopole charge of the neutron [369]. 

Although the neutron has a net zero electric charge, it is composed of charged quarks which possess a 
nontrivial radial charge distribution. This distribution produces a nonzero value of the neutron mean-square 
charge radius (r„ 2 ). To first order in the electromagnetic coupling a, the neutron-electron scattering length 

2am n c dG eN 

is proportional to the slope of the electric form factor of the neutron, G c n, in the q 2 —> limit, where q is the 
momentum transfer. For the proton (neutron), this limiting value for the electric form factor is normalized to 
one (zero). This leads to G e Ar(— g 2 ) — > ^(r 2 )q 2 , where — q 2 = q 2 is the four-momentum transfer. Although 
defined for arbitrary q 2 , in the Breit frame G c n has an interpretation as the spatial Fourier transform of the 
charge distribution of the neutron [370]. 

The sign of this slope, or equivalently the sign of the charge radius, has physical significance. For the 
neutron one expects a negative charge radius from its virtual pion cloud [371] . From the QCD point of view 
the neutron charge radius is especially interesting since it is more sensitive to sea quark contributions than 
the proton charge radius, which has a large valence quark contribution. With the advent of improved lattice 
gauge theory calculations of nucleon properties [372] and the potential to use chiral extrapolation procedures 
to ensure proper treatment of the nonanalytic chiral corrections [373] , it is possible that the neutron-electron 
scattering length may be calculable in the near future directly from QCD. Given the recent improvements 
in the measurement of the neutron electric form factor at larger momentum transfers, it is becoming more 
important to obtain a reliable value for a ne . 

There are two clusters of values in a ne measurements. One set comes from measurements of the asymmet- 
ric angular distribution of neutron scattering in noble gases, a ne = (—1.33 ± 0.03) x 10~ 3 fm [374,375], and 
the total cross section of lead, a ne = (—1.33 ±0.03) x 10~ 3 fm [376]. The other set comes from measurements 
of the total cross section of bismuth, a ne = (—1.55 ± 0.11) x 10~ 3 fm [377], and neutron diffraction from 
a tungsten single crystal, a ne = (— 1.60 ± 0.05) x 10 _3 fm [378]. Although different contributions to these 
measurements from neutron resonances and interference corrections to the total cross section from coherent 
clastic scattering have been discussed as possible sources for this discrepancy [379-381], there seems to be no 
broadly accepted consensus on the explanation for the difference. Two new experiments are in preparation; 
one exploits dynamical diffraction in a perfect silicon crystal [382], and a second attempts to improve on 
the technique of scattering in noble gases [383]. An experiment using Bragg reflections in perfect silicon 
crystals to determine a ne has also been proposed [384] . A calculation for a ne has been performed in a recent 
paper [385]. 

The precision of the charge radii of the proton and the deuteron has greatly improved over the last decade 
from theoretical and experimental advances in electron scattering and atomic physics. The charge radius 
of the proton is well-determined from both electron scattering data, Jr 2 — (0.895 ± 0.018) fm [386], and 



from high precision atomic spectroscopy in hydrogen, yV 2 = (0.890 ± 0.014) fm [387]. The charge radius of 

the deuteron is also well-determined from electron scattering data \fr 2 d — (2.128 ± 0.011) fm [388], a value 
consistent with theoretical calculations of deuteron structure from the deuteron wave function and the triplet 
n-p scattering length (-y/rj = 2.131 fm [389]), and from high precision atomic spectroscopy measurements of 

the 2P-2S transition in deuterium {yfr% = (2.133 ± 0.007) fm [390]). Atomic physics measurements of the 

H-D isotope shift of the 15-25 two-photon resonance were used to derive an accurate value for the difference 
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between the mean-square charge radii of the deuteron and proton of r\ — r 2 = (3.8212 ± 0.0015) fm 2 [391]. 
Since the neutron charge radius is simply related to the proton and deuteron charge radii, it is very timely 
for a theoretical analysis that uses these precise values as input and predicts the neutron mean-square charge 
radius in an EFT analysis. 

5.4 Neutron Polarizability 

The electric and magnetic polarizabilities of the neutron characterize how easily the neutron deforms under 
external electromagnetic fields. The quarks in the neutron are confined by the strong interaction with a 
tension equivalent to about one ton over their distance of separation of one fermi, so the polarizabilities 
are very small. To measure neutron polarizability, one may exploit the electric fields accessible on the 
surface of heavy nuclei which polarize the neutron to give a calculable contribution to the neutron-nucleus 
scattering length with a linear dependence on the neutron momentum k. A measurement using 208 Pb 
observed a term whose size and neutron energy dependence was consistent with a nonzero polarizability 
of a n = [12.0 ± 1.5(stat) ± 2.0(sys)] x 10~ 4 fm 3 [392]. Subsequent analyses have asserted that the data 
analysis is not definitive [393-395]. Another approach using deuteron Compton scattering gave a„ = [8.8 ± 
2.4(stat+sys) ± 3.0(theo)] x 10~ 4 fm 3 [396] while quasi-free Compton scattering from the deuteron gave 
a„ = [12.5 ± 1.8(stat)^ ( ! ) 'g(sys) ± l.l(thco)] x 10~ 4 fm 3 [397]. The theoretical uncertainties come from 
different treatments of strong interaction effects. 

QCD effective field theory is developing into a quantitative theory for the calculation of low energy nuclcon 
properties. The lowest-order prediction of chiral perturbation theory for the neutron polarizability is a n — 
12.2 x 10~ 4 fm 3 [398], and analysis of the extensive new Compton scattering data on the proton and deuteron 
using an EFT approach has extracted a neutron polarizability of a n = [13.0 ± 1.9(stat)^'g(sys)] x 10 _4 fm 3 
in agreement with the chiral perturbation theory prediction [399]. Lattice gauge theory calculations of 
the polarizability are also improving [400], raising the possibility of a sharp prediction for the neutron 
electric polarizability directly from QCD. This work should motivate further efforts to improve the neutron 
polarizability measurement. 

5.5 Ultrahigh Resolution Gamma Spectroscopy and Metrology with Neutrons 

A rather unexpected use of low energy neutrons in nuclear physics and absolute metrology lies in the field of 
ultrahigh resolution gamma spectroscopy [401]. The energy of a slow neutron is several orders of magnitude 
smaller than the few MeV energies of the gammas that are usually emitted upon neutron capture in a 
nucleus. As the excited nucleus emits a cascade of gammas the nucleus recoils and the subsequent gammas 
possess Doppler-shifted energies. The narrow lineshapes of these gammas, which can be resolved with 10~ 6 
precision using diffraction from perfect crystals, are sensitive to the lifetimes of the nuclear states in the 
cascade. This technique of Gamma Ray Induced Doppler Broadening (GRID) [402,403] has been developed 
at the ILL and used to identify interesting two-phonon states in various deformed rare-earth nuclei starting 
with 168 Er [404] which are of interest for understanding the degree of collectivity in these highly-deformed 
states of matter. 

The very high resolution of the crystal diffraction gamma spectrometers coupled with the special knowl- 
edge of the lattice constants of the Si and Ge crystals used, which can be calibrated absolutely in terms of the 
meter, also make possible absolute measurements of nuclear binding energies of certain nuclei by measuring 
the energies of all the gammas released in the cascade. These measurements have been performed using the 
G AMS4 crystal spectrometer at the ILL [405] . The measurement of the very small diffraction angles of the 
gammas are performed with Michclson angle interferometers of nanoradian precision. This device has been 
used to determine the mass of the neutron to a precision of 2 x 10~ 7 in atomic mass units. [406] Recent 
measurements of the binding energies of 29 Si and 33 S are in agreement with independent determinations 
using the cyclotron frequencies of ions in a Penning trap [407] and have been used to achieve the most 
precise direct test of the formula E = mc 2 , which is verified at the 4 x 10~ 7 level [408]. 

Another unexpected application of precision neutron measurements lies in the determination of the 
electromagnetic fine structure constant a. By simultaneously measuring both the wavelength and the ve- 
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locity of a neutron beam, one can measure the quantity = v\. Comparing with the Rydberg constant 
-Roo = 1 /2^ca 2 and expressing ^ = and noting that the mass ratios are known to very high 

accuracy from other measurements, one sees that a measurement of determines a. This ratio was mea- 
sured by selecting a wavelength using backscattering from a perfect silicon crystal and using polarized neu- 
trons and a resonant neutron spin flipper to measure the velocity [409]. The result of a' 1 — 137.0360113(52) 
is in agreement with and of comparable accuracy to other measurements. 

6 NEUTRONS IN ASTROPHYSICS AND GRAVITY 

This section discusses some of the ways in which neutrons and nuclear reactions involving neutrons play 
vital roles in several astrophysical processes. Neutrons play a decisive role in determining the element 
distribution in the universe. The decay rate of the neutron determines the amount of primordial 4 Hc in 
BBN theory, and neutron reactions in stars form most heavy nuclei beyond iron. In addition, one can use 
the fact that neutrons in the gravitational field of the Earth see potential differences comparable to those 
from the strong and electromagnetic interactions as an opportunity to search for gravitational effects on an 
elementary particle. 

6.1 Big Bang Nucleosynthesis 

Neutron decay influences the dynamics of Big Bang Nucleosynthesis (BBN) through both the size of the 
weak interaction couplings qa and gy and the lifetime. The couplings determine when weak interaction 
rates fall sufficiently below the Hubble expansion rate to cause neutrons and protons to fall out of chemical 
equilibrium, which occurs on the scale of a few seconds, and thus the n/p ratio decreases as the neutrons 
decay. The lifetime determines the fraction of neutrons available as the universe cools, most of which end 
up in 4 He [410], and occurs on the scale of a few minutes. The neutron lifetime remains the most uncertain 
nuclear parameter in cosmological models that predict the cosmic 4 He abundance [411,412]. With the recent 
high-precision determination of the cosmic baryon density reported by the WMAP measurement of the 
microwave background [413], there is a growing tension between the BBN prediction for the He abundance, 
which is quite sensitive to the neutrino sector of the SM, and that inferred from observation [414]. Recall 
that BBN calculations predicted that the number of light neutrinos that couple to the Z was about 3 before 
the LEP measurements ended all doubt. The quantitative success of BBN is now routinely used to constrain 
various aspects of physics beyond the SM. The small size of the baryon density relative to the density required 
to close the universe is one of the observational cornerstones of the dark matter problem in astrophysics. 

The main concern in deviations between BBN theory and experiment remains the astronomical deter- 
minations of the 4 He abundance, whose systematic errors are perhaps not yet fully understood. As the 
neutron lifetime measurements improve, other neutron-induced reactions in the early universe, such as the 
n + p — > d + 7 cross section, will become more important to measure precisely. Some reactions that are 
difficult to measure would benefit from the application of EFT methods for calculation in the relevant energy 
regime, and again low energy neutron measurements will be useful to fix the EFT parameters. It is therefore 
likely that neutron measurements will continue to be relevant for BBN. 

6.2 Stellar Astrophysics 

Nuclear reactions involving neutrons play vital roles in many astrophysical scenarios. For example, virtually 
all elements heavier than iron were made in environments inside stars and supernovae or more exotic en- 
vironments where neutron interactions dominate the nucleosynthesis. Neutrons dominate because charged 
particle rates are highly suppressed by large Coulomb barriers and because astrophysical conditions favor 
the release of large fluxes of neutrons. In addition, many of the most interesting new astronomical results, 
from the latest generation of observatories to measurements of isotopic anomalies in meteorites, are provid- 
ing views of the results of this nucleosynthesis with unprecedented detail and precision. Furthermore, more 
realistic models of astrophysical environments made possible by recent advances in computing are providing 
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new insights into the inner workings of these objects, as well as contributing to related topics such as galactic 
chemical evolution and the formation and age of our solar system. 

An examination of the observed elemental abundances in the solar system, together with rudimentary 
nuclear physics considerations, reveals that neutron capture reactions are essential for the origin of the 
elements heavier than iron [415]. Almost all these elements are thought to have been synthesized inside 
stars, supernovae, or other more exotic environments through sequences of neutron capture reactions and beta 
decays during the so-called slow neutron capture ("s") [416] and rapid neutron capture ("r") [417] processes. 
The s and r processes are each responsible for roughly half of the observed heavy element abundances. 
The remaining neutron-deficient isotopes that cannot be reached via neutron capture pathways are thought 
to have been formed in massive stars or during supernova explosions through the photodissociation ("p") 
process. 

In some cases, further progress in these areas is hampered by the lack of accurate rates for nuclear 
reactions governing stellar nucleosynthesis. Many of these astrophysical reaction rates can be determined by 
measuring neutron-induced cross sections in the energy range between approximately leV and 300 kcV. For 
about 20 radionuclides along the s-process path, the neutron-capture and /3-decay time scales are roughly 
equal. The competition between neutron capture and (3 decay occurring at these isotopes causes branches in 
the s-process reaction path that, if measured, could be used to directly constrain dynamical parameters of 
s-process models. There is very little data on the (n,7) reaction rates for such radioactive branching points. 

Measurements of cross sections in an energy range relevant to astrophysics have been conducted at 
facilities such as ORELA [418], GELINA [419], nTOF [420], DANCE [421], and IREN [422] along with 
others and could also be pursued at new epithermal neutron sources such as the SNS and JSNS [423] . The 
intensities now suffice to conduct cross section measurements on small quantities of unstable isotopes. As the 
understanding of the origins and processes which lead to clement formation improves, we can more effectively 
exploit astrophysical observations to constrain the understanding of what phenomena may lie beyond the 
SM. 

6.3 Gravitationally-induced Phase Shift 

The contribution of precision neutron measurements to gravitation are few in number but notable in con- 
ceptual impact. The equivalence principle for free neutrons has been verified at the 10 -4 level [424], and 
although one might do better with ultracold neutrons [425], experiments with bulk matter are several orders 
of magnitude more precise. Another connection between neutron physics and gravity is the observation of 
the gravitational phase shift by neutron interferometry, which was the first verification that the principles 
of quantum mechanics seem to apply to the gravitational potential as well as the potentials produced by 
other interactions [426,427]. This measurement has been performed with an accuracy at the 1% level, and 
there is a slight disagreement between theory and experiment [428] . The source of the difference is believed 
to lie in gravitationally-induced distortions in the perfect crystal interferometer as it is rotated to change 
the relative height of the two paths in the interferometer. There are two plans underway to improve the 
measurement, one to conduct the experiment with the interferometer suspended in a neutrally-buoyant fluid 
to eliminate possible distortions of the interferometer crystal [429] and another to use a recently-developed 
Mach-Zehnder interferometer with cold neutrons [430]. 

6.4 UCN Gravitational Bound States 

Experimental tests of gravity and searches for new long-range forces have attracted more interest in recent 
decades. A reanalysis [431] of an experiment on the principle of equivalence by Eotvos motivated a series of 
precise tests of the principle of equivalence culminating in the torsion balance experiments of the Eot-Wash 
group, which set stringent new limits to violations of the equivalence principle [432]. Recently, speculations 
involving the propagation of the gravitational field into extra dimensions produce as a natural consequence 
a modification of the inverse square law for gravity on submillimeter scales [433] . The current experimental 
limits on the existence of such forces are not stringent enough to rule out this possibility. Experiments that 
use smaller separations in torsion-balances [434] and small mechanical resonators [435] have already set useful 
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limits. Such experiments share some techniques with those used in the recent experimental demonstration 
of the Casimir effect [436]. 

At first glance experiments with neutrons would not seem to offer a productive technique to conduct 
sensitive searches for new short-range forces of gravitational strength. However, the very small polarizability 
of the neutron is an advantage relative to atoms, whose van der Waals attraction poses a background issue 
for such searches. Although neutron beam densities are very small compared to bulk matter, the delicate 
control possible for neutron polarization and the weak interaction of the neutron with material media make 
it possible to pass the neutron through bulk samples, thus presenting an interesting opportunity to probe 
new forces at a distance scale that is difficult to reach with other methods. Possible experiments to probe 
extra dimensional gravity theories using the angular and neutron energy dependence of neutron scattering 
from spin zero nuclei have also been discussed [437]. 

The recent measurements in search of gravitational bound states of ultracold neutrons are of some interest 
for constraints on new forces at smaller distance scales [438,439]. For neutron kinetic energies small compared 
to the neutron optical potential of a plane surface, the potential seen by a neutron moving above the plane 
in the gravitational field of the Earth should possess neutron bound states with energies given to good 
approximation by the zeroes of Airy functions, which are solutions to the Schrodinger equation in a linear 
potential. The lowest bound state of energy 1.4 peV hovers above the medium at a distance of order 10 /um, 
and any nonstandard attractive interaction of the neutron with the matter on this length scale could create 
another bound state. 

The experiments were designed to populate these bound states by forcing the UCN to pass through a 
narrow gap above a planar medium and to detect their presence by measuring the transmission of the UCN 
through the gap as the separation is varied. Intensity oscillations in the transmission are observed which 
can be fit to a model which includes the effect of the spatial extent of the bound states. The agreement of 
this data with theory was used to set an interesting bound on gravitational-strength forces in the nanometer 
range [440,441]. 

Other sources of constraints on new forces in this distance regime from neutron measurements come from 
comparisons of the results of different experiments to measure neutron-electron scattering lengths. Total 
cross section measurements are sensitive to the forward scattering amplitude through the optical theorem, 
but measurements of the angular distribution of the differential cross section depend on the momentum 
transfer. In the presence of a new force on the distance scale set by the momentum transfer, the neutron- 
electron scattering amplitude inferred from the two measurements could be different. The agreement of the 
neutron-electron scattering length from these two different measurements sets a competitive limit on new 
forces at the nanometer scale. Analysis of total neutron cross section measurements on 208 Pb using eV to 
keV energy neutrons can also be used to set a limit for shorter distances [442] . 

7 SUMMARY 

Slow neutron experiments address fundamental scientific issues in a surprisingly large range of physics sub- 
fields. The development of new types of ultracold neutron sources, pulsed and CW spallation sources, and 
the continued increase in the fluence of reactor beams is making possible new types of experiments and 
opening new scientific areas. We anticipate significant progress in the areas of neutron decay measurements 
and neutron EDM searches and also in the field of weak NN interactions and the measurement other low 
energy neutron properties. We also expect more use of neutron measurements for applications in astro- 
physics. As is typical, many of the new opportunities made possible by technical developments and the 
emerging scientific issues were not clearly foreseen a decade ago. The steady progress over the last decade 
in quantitative theoretical understanding of the strong interaction is an underappreciated development with 
exciting applications to neutron physics. 

Although still a relatively small field compared to other areas in nuclear and particle physics, the ex- 
panding scientific opportunities in fundamental neutron physics are attracting a growing number of young 
researchers. This growth is driving the increasing number and variety of new facilities. The diverse appli- 
cations and the location of neutron physics at facilities whose main purpose is generally neutron scattering 
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have combined to obscure somewhat its accomplishments. We hope that the reader of this review has gained 
an appreciation for the breadth of activity in the field. 
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